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The  Distribution  of  Space  Charge  in  the  Hull  Magnetron  Diode 

by 

John  A.  Bradshaw 

Cruft  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 

Abstract 


If.  in  a magnetron  diode,  the  current  is  not  limited  by  cathode 
temperature,  then  the  radial  dependence  of  the  potential  and  the  charge 
density  functions  may  be  found  before  cutoff  from  Langmuir's  solution, 
modified  by  effects  of  increasing  transit  time.  In  the  cutoff  transition 
and  after  cutoff,  the  cathode  temperature,  as  reflected  in  the  distribution 
of  electronic  velocities  at  emission,  affects  the  potential  and  density  func- 
tions. Data  on  anode  current  characteristics  are  compared  with  cal- 
culations based  on  models  tliat  include  these  temperature  effects. 

The  data  from  T.R.  135  on  probe  signal  absorption  are  interpreted 
as  showing  that  the  space-charge  cloud  extends  clear  to  the  anode  even 
after  cutoff,  and  is  not  bounded  by  the  so-called  Hull  radius  . This  and 
other  evidence  point  to  electronic  interactions  as  the  source  of  a 
mechanism  for  traneporting  charge  through  the  cloud  from  cathode  to 
anode.  Finally,  the  electronic  distribution  in  phase  space  is  stated  in 
a form  which  could  be  verified  and  refined  by  further  experiments. 


The  Distribution  of  Space  Charge  in  the  Hull  Magnetron  Diode 

I 

Transit  Time  Effects  before  Cutoff 

a.  Introduction 

This  report  supplements  Technical  Report  No,  185,  entitled  "A 
Probe  Signal  Study  of  the  Hull  Magnetron  Diode."  In  the  following  pages 
we  attempt  to  fit  the  TEM  probe  resonance  data  into  a self-consistent 
model  of  the  magnetron  diode.  In  the  fifth  and  last  chapter  we  present 
the  most  controversial  aspects  of  our  model  of  the  space  charge  distri- 
bution after  cutoff.  In  the  chapters  preceding  that  one  we  present  further 
data  on  emission  and  leakage  currents  in  the  experimental  diode,  and  ex- 
amine in  some  detail  the  effects  of  a thermal  distribution  of  emission 
velocities  on  the  anode  cui'rent  and  on  the  space  cliarge  distribution.  In 
this  first  chapter  we  consider  the  effects  of  the  transit  time  cf  electrons 
in  a weak  magnetic  field. 

The  Hull  cutoff  transition  in  the  plate  current  characteristic  of  a 
magnetron  diode  clearly  divides  the  problem  of  the  radial  distribution  of 
space  charge  within  the  tube  into  two  sections.  Since  the  main  features 
of  these  sections  and  the  boxxndary  between  them  depend  chiefly  on  the 
ratio  p^of  steady  anode  potential  to  the  square  of  the  rragrieiic  field 
imposed  on  the  tube,  we  may  follow  the  French  usage  and  call  these  two 
sections  regimes,  the  regime  before  cutoff,  and  the  regime  after  cutoff. 
We  thus  reserve  the  word  region  for  spacial  regions  within  the  diode. 


The  regime  before  cutoff,  which  we  consider  in  this  chapter,  is 
fairly  well  understood.  At  the  extreme  p — >oo,  we  have  the  Langmuir 
solution  for  a diode  without  magnetic  field.  If  a magnetic  field  is  im- 
posed, and  if  p decreases  to  the  critical  or  cutoff  value,  p^,  the  increase 
in  the  transit  time  of  individual  electrons  serves  to  explain  the  observed 
decrease  in  diode  current.  The  effects  of  transit  time  were  calculated 
by  Benham  and  Bethenod  tor  a planar  magnetron,  and  by  Allis  andBrillouin 
for  a cylindrical  one.  We  develop  an  integral  e<^uation  as  an  alternative 


1 


and  perhaps  more  convenient  way  to  find  the  anode  current  as  function  of 
p in  a cylindrical  tube. 

b.  Two  Calcvilations  for  the  Linear  Magnetron 


In  a diode  without  magnetic  field  the  thermionic  electrons  move  to  the 

anode  along  practically  straight  lines.  Their  paths,  it  seems,  depend  little 

on  their  emission  velocities,  and  their  number,  outside  the  virtual  cathode, 

depends  little  if  at  all  on  the  total  emission.  Imposition  of  a weak  magnetic 

field,  of  course,  deflects  these  electrons  so  that  they  trace  longer  curved 

paths,  spend  more  time  in  transit,  and  hit  the  anode  at  an  oblique  angle.  The 

anode  current  falls  slightly,  and  the  reason  is  briefly  this:  the  virtual 

cathode  supplies  just  enough  current  to  maintain  a complete  shield  (£^=0) 

against  the  anode  potential.  The  longer  an  electron  spends  in  transit,  the 

longer  its  charge  helps  to  shield  the  cathode.  Hence  if  the  average  transit 

time  Tj  increases,  fewer  electrons  need  to  be  released  per  unit  time  to  main- 
d 

tain  the  shield,  and  consequently  falls. 

For  estimating  the  transit  time  and  its  effects,  W.  Benham  [1]  in  1935 
published  an  extension  of  the  Langmuir 'Child's  equation,  which  he  had  worked 
out  for  a planar  magnetron.  Following  him,  we  assume  a stead/  state,  so 
the  total  current  per  unit  area  J^,  across  any  plane x parallel  to  the  planes 
of  the  cathode  and  anode,  x = 0 and  x = d respectively,  is  independent  of  time 
and  of  X.  The  electrons  have  actually  a distribution  in  velocities  given  by 
some  density  function  f(x,x),  so  we  write  and  the  corresponding  charge 
density  p(x)  as: 


CD 


J = e 


/*  iU(x,x)dx 


-GO 


p(x)  = e J 


f(x,x)dx 


X = 


dx 

•at 


(1.1) 


-oo 


For  Ej  the  electric  field  within  the  diode,  we  have  the  following  rela- 
tions, the  first  by  assumption,  the  second  as  Poisson's  equation,  and  the 
third  Eviler's,  with  the  path  parameter  <T  ; 


8E  _ 0 ^ _ P(x)  ^ _ 8E  ^ 

■ ’ 8x  € ’ dT*  “ "Bx  d'T 


(1.2)  . 


Defining  a velocity  x,  an  average  over  the  electronic  distribution,  by 
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= p(x)  x(x),  we  conclude: 


^ - dx  • J.I.  ciE  ‘ X . . 

If  X , then  ~ t a constant. 


(1.3) 


With  equation  (1.3),  often  called  Llewellyn's  equation,  one  may  integrate 
the  equations  of  the  motion  of  an  average  electron  in  a planar  diode.  How- 
ever, this  reqvires  the  average  electron  to  be  typical  too,  or  the  deviations 
from  the  average  to  be  negligible.  For  interpenetrating  beams  of  particles, 
or  even  for  swarms  with  a wide“«pread  in  random  velocities  such  as  make  tip 
a virtual  cathode,  Llewellyn's  equation  has  important  limitations.  These 
restrictions  are  discussed  by  J.  R.  Knipp  [2,  page  135],  and  apply  partioilarly 
to  the  magnetron  after  cutoff,  (Chapter IV). 

To  continue  Benham's  method  we  assume  a uniform  magnetic  field  B, 
directed  along  the  z-axis,  and  represented  by  ci)  = eEAn.  The  average  elec- 
tron, released  from  the  virtual  cathode  at  x = 0 with  zero  velocity  at  time  t^, 
obeys  the  force  equation; 


^ 2 eE  _ eJ  _ d x 

x"  + X = —=— -I  = C7^  , X ' = — ,r 
m me  ^ ^ 


(1.4) 


This  relation  implies  E(x=0)  = 0.  Two  integrals  of  it,  with  the  initial  con- 
ditions x'(0)  = x"(0)  = 0 are: 


x'  = (1  - cos  0)  , X = (0  sin  0)  . 

(Jj  U) 


(1.5) 


We  call  6 the  transit,  angle.  If  ^(x)  is  the  potential  relative  to  the  virtual 
cathode,  a general  first  integral  of  (1.4)  is  the  energy  equation: 


(1.6) 


Substituting  from  (1.5)  in  (1.6)  for  x'  and  x,  evaluated  at  x = d,  we 
find  p,  the  normalized  potential  at  the  anode,  as  a function  of  the  angle  of 
transit  to  the  anode,  6^.  Giving  the  current  density  J a similar  normaliza- 
tion j,  we  have: 


1 + j^(l  - cos  0^)^j  , j = (0^  - sin  e^)"^ 


(1.7) 


■p  T > 
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Here  by  definition, 


P = - 


eV 

a 

2TI 

mu  d 


J = 


eJ 

— 

mcu  d 


With  (1. 7),  one  could  plot  j(p)  with  0^  as  a parameter.  (Brillouin  does 
this  [3,  page  95]  without  a reference  to  Benham. ) In  Fig.  1-1  we  plot  in- 
stead the  ratio  R{p): 


(planar)  "TT 


9 j 


4S'Tp 


37Z- 


(1.8) 


Rp  shows  how  much  a small  magnetic  field  reduces  the  anode  current  1^  in 
a planar  diode  from  its  Langmuir  value  In  the  limit  0^  = 0,  or  p->co 

we  have  the  Langmuir -Childs'  law,  R = 1,  In  the  limit  6 = or  p = p^,  the 
electronic  orbits  just  graze  the  anode  and  experimental  curves  of  the  cutoff 
transition  begin  to  show  the  effects  of  a thermal  spread  in  emission  velocities 
The  value  R(p^.)  = 0.  717  (Brillouin,  page  98  of  [3])  has  little  significance.  In 
the  regime  well  before  cutoff,  however,  where  Bcnham’s  assumptions  apply, 
his  method  gives  ua  a complete  picture  of  the  planar  magnetron  and  of  the 
effects  of  on  the  anode  current. 

One  may  calculate  a relation  between  p and  j without  introducing 
Llewellyn's  equation  explicitly.  Bethenod  in  1939  [4]  gave  a method  which 
we  sketch  here  briefly  because  it  introduces  the  method  we  use  for  a 
cylindrical  magnetron,  although  his  results  are  less  transparent  than 
Benham's.  He  takes  as  basic  relations  the  energy  equation  (1. 6),  Poisson's 
equation,  and  the  current  relation: 


d _ p 


= X p(x)  = constant 


(1.9) 


He  regards  x rather  than  4* as  the  independent  variable. 
We  introduce  again  normalizing  definitions: 

e J 


s = (x,'d)  , ^(s)  = 

mw  d 


me  dw 


Then  using  (1.6),  (1.9)  and  an  integration  by  parts  we  have; 


I 
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i(2#(s)  - 8 


s 

*(»)  = j J 


(s-a*)  da* 


(1.  10) 


8*)  - (8*)‘ 


In  the  formal  integration,  8 * is  the  integration  variable.  The  integral  form 
is,  for  the  present,  a blind  alley. 

2 2 

Bethenod  introduces  an  effective  potential,  "y  = ^ -f  — — , a formula 
which  Twisa  [5]  and  many  others  have  used,  Some  writers  interpret ')^(x) 
as  the  potential  in  a coordinate  system  moving  at  the  velocity  y(x)  = sue  in 
the  direction  perpendicular  to  both  the  B-  and  the  E-fieldi. . The  value  of 
this  interpretation,  and  the  convenience  of  the  formiila  is  often  open  to  ques- 
tion, but  in  Bethenod' 8 problem  it  allows  us  to  complete  the  integration. 

Accordingly,  we  define  f(s)  to  the  fourth  power,  to  avoid  f.ractional  ex- 
ponents: 


f^(s)  = 24>(s)  - 8^ 


(1.11) 


With  inititt).  conditions  f = 0 = df^/35  at  s = 0,  the  integrals  of  (V  \l) 


- 8jf^  - 4f^ 


f(s)  2 

. = 2 r 

o 


(1.12) 


Writing  for  convenience,  y = f(2j)”  , we  find  ai  s = 1: 

"d  2 , , , 

1 = 4j  / -y  = 2j  [sin'Vd  - " ^d^^ 

^o  All 


(1.13) 


1/4 

This  is  an  implicit  relation  between  p and  j,  since  we  have  y,  = (2p-l) 

J2j)”  The  same  expression  (1.  13)  can  be  obtained  from  Benham's  equa- 
tions (1.7)  after  eliminating  between  them  the  angle  0^. 

c.  An  Integral  Equation  for  the  Cylindrical  Magnetron. 

In  the  cylindrical  diode,  Llewellyn's  equation  does  not  yield  integrable 
relations  for  the  electron's  average  orbit.  To  show  this  we  give  the  radial 
equations  comparable  with  (1.3)  and  (1.4): 
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drE 

'dr 


2v<  n 


■—  - V 


(1.14) 


Here  h Is  the  axial  length  of  the  cathode.  The  uniform  magnetic  field  B 
lies  parallel  to  the  cathode  and  to  the  s-axis,  and  is  represented  aghin  by 
u = (eBi^).  E(r)  is  the  radial  electric  field,  derivable  from  a potential 
function  d(r),  1^  is  the  total  radial  current,  a constant,  and  P,  also  a cisa-v 
Btant,  is  the  normalized  form  of  the  canonical  angtilar  momentum.  If  the 
radius  of  the  cathode  is  r and  we  introduce  the  normalization  s = (r/r  ) , 

C V 

then: 

P=i£(0+“)  , (TR185.  1.6)  (1.15) 

w 2 

ft  a 

If  the  emission  velocities,  r and  r8  at  r = r^,  are  negligible.  MM  we  as.* 
sumed  above  for  the  linear  magnetron,  then  P = 1. 

Allis  in  1941  [6,  pages  10  40]  developed  a solution  of  (1.  14)  by 

inserting  in  it  a series  (r/r^.)  = ^ and  evaluating  the  a^  coefficients. 

o 

Slater  comments  [?,  page  3 53]  that  these  calculations  of  Allis  "do  not 
seem  to  be  particularly  accurate"  in  the  regime  before  cutaj!!f.  The  seven 
points  on  the  curve  R_(p)  = 1_/1t  » curve  for  the  cylindrical  diode  in 
Fig.  1-1,  were  therefore  calculated  by  anew  method.  We  devote  the  reft 
of  this  chapter  to  an  outline  of  it. 

The  method  is  based  on  iteration  of  an  integral  equation  similar  to 
the  integral  in  (1.  10).  We  assume  a potential#^,  the  normalised I<angiavir 
potential  which  is  correct  in  a diode  without  magnetic  field.  We  calculate 
in  such  a potential  the  transit  time  7“ as  modified  by  a small  magnetic 
field.  We  then  calculate  the  potential  perturbed  by  a modification  of  the 
transit  time.  These  two  steps  will  be  somewhat  obscured  by  an  integra- 
tion by  parts.  They  could  be  repeated  to  improve  the  potential  further. 

One  could  start,  we  note  parenthetically,  with  a charge  free  diode 
in  a magnetic  field,  and  the  appropriate  (logarithmic)  potential.  One  would 
then  treat  a small  current  in  this  diode  as  a perturbing  influence  on  the 
potential.  Alter  several  iterations  one  should  arrive  at  the  same  reexilt; 
the  potential,  the  current  and  the  transit  time  adjusted  into  one  aelf-consiatent 


3 

J 

i 

J 


;■  ■ 
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solution.  However , is  closer  than  a logarithmic  form  to  the  desired 

potential  ¥(s);  an  integration  over  should  therefore  give  a closer  first- 
order  result. 


To  find  the  equation  for  a cylindrical  diode  that  will  correspond  to 
(1.  10)  for  a planar  one,  we  take  again,  this  time  in  cylindrical  coordinates, 
Poisson's  equation,  the  energy  equation,  and  the  current  relation; 


dr  d^  __  ^ 

r 3F  3F  ~ c 


r"  + (r6) 


m 


= Zirhrrp 


(1.16) 


We  assume  the  virtual  cathode  at  r = r^,  the  reference  for  the  potential 
releases  electrons  with  negligible  velocities.  Then  P is  indeed  one,  and 
{1.  15)  may  be  used  to  eliminate  6.  We  rewrite  (1.  16)  in  terms  of  a and  a 
normalized  potential^: 


d sd4»  . 

cTi'^T'  ^ 


4mw  « 


s^  + w^{s-l)^  = 8u^s® 


I = irh  s pr 
a c 


(1.17) 


We  find  the  following  definitions  convenient  throughout  this  report: 

reE(r)  el^ 


mw  r 

c 


•; — r~T 

2mu>  r 

c 


M = 


TZ — 5 T 

4irhmu  c r 
o c 


(1.  18) 


Here  K we  regard  as  the  normalized  electric  field  zind  M £.s  the  normalized 
radial  current,  a constant.  If  the  electric  field  at  the  cathode  vanishes, 

K(r  ) = 0,  and  the  equations  like  (J.  105 


V • 


I 


d d4> 
ds  * ds 


M 


^8a9-(B-ir 


or  <I>(b 


' - M r (1^  s - In  s* ) ds* 
n ‘^8s*«>(s»)-(s*-l)^ 


(1.19) 


Since  we  propose  to  build  our  iterated  solution  fr  om  the  >7ell  known 
Langmuir  solution,  v.'e  must  write  Langmuir's  solution  in  these  symbols. 

If  tie  magnetic  field  vanishes,  <S>  prove  a inconvenient,  so  for  two  paragraphs 
we  drop  this  normalization  and  write  »)<  = " e^/mr^^  = f>.  Then,  after 

multiplying  the  first  equation  in  (1*19)  by  4s,  we  find: 

s£e£ — •„  5 = >in, 


.2  , 

a _ 


2 ‘*rh  r <m 
c 


•5? 


(S'l)' 


(1.20) 
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In  (1. 20)  we  may  set  w = 0.  Then,  writing  (L20)into 

two  relations,  one  which  normalizes  and  one  which  givp^s  as  a func- 
tion of 

^ T”  . (A  > ^ • (1-21) 

8 tr  «h  “m 

Langmuir  and  Motl-Smith  [3]  developed  fj^  in  a MacLiaurin  series  in 
§from  the  equation  in  (1. 21)  for  f(^),  and  from  derivatives  of  this  equa- 
tion, all  evaluated  at  the  cathode  where  % = 0.  For  instance  4<(0)  = 0 =f{0), 
and  f^(0)  = 1 from  {1.  21).  We  find  from  the  first  derivative  of  (1.  21): 

+ 3fj^f"j^  = expCg)  , fj^(0)=i  . (1.22) 


This  process  continues  to  yield  one  coefficient  per  derivative,  and  the 
resulting  aeries  is; 


•••  = 


exp  (5/2)  • p(^) 


(1.23) 


fl(5)  is  the  function  usually  quote.*!  as  Langmuir's,  rather  than  and 
Langmuir  showed  that  p approaches  one  as  ^ tends  to  infinity.  The  series 
for  p is  [9] : 

P h’§  - I . ..  (1.24) 

From  (1.  20)  with  = 0,  we  conclude  that  the  equations  which  f^ 

Xi 

satisfies  can  be  written  in  terms  of  s as  follows; 


d 

Ts'  “ 


df. 


4/3 


da 


9Wsf 


or 


£l(s) 


4/3 


S 

1 r (In  - “ In  8*  \fi«t  4c 


'1  A/s*f 


(i.  25) 


(s*) 


From  (1.  23)  we  concluUe  that  the  solution  of  these  is: 


r 


- T-  - 
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magnetic  field,  by  comparison  of  (1.  25)  with  (1,  19): 

s 

,4/3,  v_l  r (In  8 -In  8*)  ds 

r 4/y ' — ’ 

y8*f^^^8*)-S{8^,8*) 


(1.27) 


When  we  have  found  the  function  S in  (1. 27)  and  brought  the  current  density, 
3/2 

M in  (1.  19)  or  k in  (1. 21),  into  the  iterated  solution,  the  formulation  will 
be  complete . 

Now  the  expression  imder  the  square  root  in  the  denominator  of  each 
integral,  (1.19).  (1.20),  (1.25)  or  (1.27),  represents  the  radial  velocity 
s ^(1.  17).  In  T.R,  185  (1.  5)  we  defined  for  orbits  in  a strong  magnetic 
field  a turning  radius  s^^,  where  s vanishes.  In  the  regime  before  cutoff , the 
average  electron  hits  the  anode  before  reversing  its  radial  velocity,  but  we 
may  imagine  its  path,  and  the  potential  f,  (s),  extended  beyond  the  anode 
through  the  interval  8^<l8  <Sj^.  Now  S represents  the  effects  of  a small 
magnetic  field  on  s,  so  from  (1. 17)  S must  be  proportional  to  (s-1)^,  and  S 
in  (1. 27)  must  cause  the  denominator  of  the  integral  to  vanish  at  s^: 

v2 


» (8-1) 


/4 


( 1 . 28) 


-2/3 


As  a check  on  formula  (1.28),  we  note  that  S approaches  zero  like  s^ 
as  8^  approaches  infinity  with  s fixed.  As  s^  approaches  infinity,  the 
magnetic  field  tends  to  zero  and  fj_^(8|^)  behaves  like  s^^^(1.26). 

For  the  first  iteration  in  (1.  19)  we  define  the  Langmuir  trial 
potential  by  arguments  like  those  for  S’. 

*L<”>  =*b  ■ *b  = 

This  (s)  causes  the  denominator  in  (1.  19)  to  vanish  at  s^.  The  first 
iteration  from  (1.  19),  with<t(s^)  = (1’27)  yields: 


<i>(s)  = M . 9f‘^^(a)  . 


(1. 30) 


But  if<i>(8)  is  to  be  used  for  another  iteration,  it  must  have  the  same  form 
as (s).  Hence: 


.S 


! ! 

I 

I 

i 


, ! 


.1 


T 
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r-4/3. 


4.(8)  =4.J^ 


p =*(*,) 


M =4>, 


3/2  --V3, 


<*b'  'L^<‘b' 


(1.31) 


(1.32) 


With  a little  more  algebra  one  can  find  the  quantity  R^,  plotted  on 
the  right  in  Fig  . 1 • 


•1: 


T?  = 
C 


W 


9M 

2V2 


_3/2 


(1.33) 


The  ratio  {^(■•p^lj_,(  ®v,'|  1^^®  fh®  value  1.  19  > found  by  numerical  work  to 

be  independent  of  in  the  range  7 1/4  29  1/4.  The  corresponding 

transit  angle  0^(®{j)  i«  given  by  an  integral  which,  like  (1.  27),  can  be 
evaluated  numerically,  once  Sj^  and  s^  are  given: 

. L 4/3, 

"K/  / ds 

.(1.34) 


W “J  .r— — -7  “*b  ' •-( 


1 


VS 


(s)  - S(s.Sj^) 


'This  angle  appears  in  the  parameter  ma  60  along  the  curve  R in 
Fig.  3-1. 

This  completes  our  analysis  of  transit  time  effects.  We  delay  the 

presentation  of  data  on  the  ratio  I to  Chapter  III,  where  we  consider 

the  effects  of  emission  velocities  on  this  ratio.  If  one  holds  V fixed  and 

a 

increases  w,  then  as  p falls  toward  the  cutoff  transition  value  p^,  the  cir- 
cumferential component  of  the  emission  velocity  begins  to  affect  the  path 
and  transit  time  of  the  individual  particle.  In  the  cutoff  transition  the 
sign  of  this  component  practically  determines  whether  the  particle  reaches 
the  anode  or  returns  to  the  cathode.  For  p somewhat  larger  than  p^,  how- 
ever, the  data  support  the  analysis  leading  to  Fig.  1-1,  and  the  conclusion 
that  transit  time  effects  are  the  chief  cause  of  the  observed  decrease  in 
I^  accompanying  an  increase  in  w. 


L. 
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Effects  of  Emission  Temperature  in  a Diode  without  Magnetic  Field 
a.  Introduction 

In  this  chapter,  by  a review  of  the  characteristics  of  a virtual 
cathode  without  magnetic  field,  and  by  an  analysis  of  experimental  data 
on  the  emission  current  I^,  we  prepare  the  way  for  Chapter  III,  where 
we  vdll  consider  the  effects  of  a magnetic  field  on  the  virtual  cathode, 
and  the  effects  of  emission  temperature  on  the  cutoff  transition  in  the 
magnetron  diode . 


TR201 


-12- 


The  prcblem  of  describing  space  charge  in  front  of  an  emitting 
surface  is  at  least  forty  years  old.  In  1914  W.  Schottky  [l]  published 
an  analysis  of  current  flowing  from  an  emitting  wire  to  a concentric 
cylindrical  anode  against  a retarding  potential.  He  found  a Maxwelliem 
distribution  of  emission  velocities.  In  another  important  paper  in  1923, 
I.  Langmuir  [2]  reviewed  a series  of  calculations  on  the  formation  of 
a virtual  cathode,  and  published  cables  applicable  to  a planar  diode. 

In  1940  E.  Wheatcroft  [3]  gave  the  results  of  a unified  group  of 
calculations  for  the  position  and  depth  of  the  potential  minimum  in  front 
of  a cylindrical  emitting  surface.  However,  hi*  vraphs  do  not  include 
results  applicable  to  a tube  with  a large  cathode.  We  have,  therefore, 
used  his  analysis  in  somewhat  esvised  form  to  extend  the  range  of  these 
results  so  as  to  apply  them  to  our  diode.  This  revision,  outlined  below, 
proves  adequate  to  treat  the  more  general  case  of  the  diode  in  a mag- 
netic field,  so  it  forms  the  basis  of  the  nevt  chapter  as  well. 


b.  Space  Charge  Inside  the  Virtual  Cathode. 


We  assume  the  conventional  polar  coordinate  system  (rdz);  the 

z-axis  lies  along  the  center  of  the  set  of  two  coaxial  conducting  cylinder  s 

that  constitute  the  cathode  and  anode  of  our  diode.  Our  reference  is  the 

emitting  surface  on  the  cathode  at  the  radius  r = r^.  The  potential 

fimction  )4(r)  is  zero  there,  and  for  r>  r_,  ^ is  determined  by  the  anode 

c 

voltage,  and  the  charge  distribution  which  we  wish  to  calculate.  We  as- 
aiiTTiA  in  this  chapter  that  electrons  suffer  neither  collisions  nor  indi- 
vidual interactions,  indeed,  that  the  orbits  after  emission  are  determined 
entirely  by  the  emission  velocities  and  by  through  the  energy  and 
momentum  relations: 


2 

mv 


2 

mv 

c 


2 


(2.1) 


We  v/ill  re-examine  these  assumptions  in  Chapter  IV. 

We  now  normalize  equations  (2. 1)  by  dividing  them  by  kT,  where 
k is  Boltzmann's  constant,  and  T is  the  cathode  temperature  in  degrees 
Kelvin.  We  also  adopt  the  definitions: 


Z , 1/2  . li 

= (— ) > ‘i*  = 5t  ’ ^ ^2kl^  ’ y - "®  ^2kf^ 


{2.2) 


Thus  X and  y are  normalized  velocities  at  some  r ^ r , while  the  normalized 
emission  velocities  we  write  without  bars.  The  coordinate  z is  suppressed, 
for  we  assume  that  ^ and  z are  independent  of  z and  of  time.  Equations  (2.  1) 


become: 


+ »j»  = x^  + ,V^  y = y • 


(2.3) 


If  the  distribution  of  velocities  among  electrons  on  emission  is 
Maxwellian,  then  at  the  cathode  the  increment  to  the  normalized  charge 
density  dS  in  any  increment  dxdy  of  the  plane  of  emission  velocities  (x,y)is: 


dS  = ^ exp  |^“(x^+y^)|  dx  dy 


(2.4) 


The  normalization  here  is  such  that  the  integral  S(r^)  over  the  whole  (x,y)- 
plane  is  equal  to  one.  If  I is  the  emission  current,  the  increment  to  the 
total  current  in  the  outward  radial  direction  is  similarly: 


dl(r  ) = — ^ 
c'  ,1 — 

Vw 


o f , 2 . 2.1  ^ , 

— exp  1 -\x  ry  )i  xdxdy 

TT  ^ ) 


(2.5) 


The  integral  of  dl  over  the  half  plane  x>  0 is  just  I^.  Under  our  assump- 
tions, each  increment  dl  is  prevented  from  exchanging  electrons  with  other 
increments,  so  as  dl  moves  away  from  the  cathode,  it  is  conserved: 

d(dl)  _ wish  next  to  determine  how  the  corresponding  charge  in- 

dr 

crement,  dS(r),  changes  as  r increases. 

To  examine  the  movement  of  charge  away  from  the  cathode,  we  change 
to  coordinates  (x,  y)  at  some  fixed  radius  s.  From  (2.  3)  we  have  the  rela- 


tions: 


d(x^)  + 2y  dy  = d('.«  ) + 2y  dy  , Vs  dy  = dy  . 

2 —2 

The  Jacobian  of  a transformation  from  (x  ,y)  to  (x  ,yj  is: 


(2.6) 
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B(x^) 

1 f{x.s) 

By 

■lx. 

0 Vs 

»(?^) 

37 

d(x^)  dy  =Vs  d(x^)dy 

Now  by  (2.  3)  and  (2.  8)  dl(r)  equals  dl(r^)  ’f  we  have; 


(2.7) 


(2.8) 


cL(r)  = 21  \f^ Bxo  l-(x^ 

o ITT  * L • 

With  a little  algebra  we  obtain  the  relation  of  dl(r)  to  dS(r): 

dl(r)  = 21^  Vbtt  xdS(r)  . (2  10) 

In  the  limit  r = r^  this  relation  satisfies  (2.  4)  and  (2.  5).  If  we  write  the 
unnormalized  increment  of  charge  density  in  dS  as  dp,  then  (2.  10)  represents 
the  relations; 

dl  = 2ffhrrdp  . dp  = dS  . (2.11) 


+V^+llll 


i|  X dx  dy 


(2.9) 


Equation  (2.  10)  implies  the  result  for  dS(r): 

dS(r)  exp  |-(T^+y^+4»)|  dx  dy  (2.12) 

The  charge  density,  normalized  with  respect  to  the  density  at  the 
cathode,  is  then  an  integral,  S(r)  = JJdS{T,xY)  over  a region  of  the(x,y)- 
plane.  We  call  this  region  W ard  proceed  to  determine  its  boundaries: 

If  an  electron  in  its  orbit  away  froni  the  cathode  reaches  a radius  r and 
thus  contributes  to  S(r),  its  radial  velocity  must  be  positive  from  r^  all 
the  radial  interval  to  r;  in  symbols,  we  require  x(r)>  0,  Secondly, 

if  its  total  velocity  is  insufficient  to  carry  it  over  the  maximum  of  v(<, 
i}j^  (which  is  the  minimum  of  ^(r),  since  e is  negative)  at  the  virtual 
cathode,  the  electron  reverses  its  radial  motioja  and  falls  back  to  the 
cathode  surface,  it  then  contributes  twice  to  S(r),  r<lr  These  two 
observations,  together  with  equations  (2.3),  serve  to  fix  the  boundaries 
of  W. 
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The  region  of  W which  ri»pre«ent«  electrons  falling  back  toward  the 

cathode  lies  properly  in  the  lower  half  of  the  (x,y) -plane.  Since  the  dis- 

r —2  —2  T — 

tribution  function  exp  •j^-(x  +y  +4,-)j  is  symmetric  in  x we  may  reflect  this 
region  about  the  y-axis  into  the  upper  half  of  the  (xy)-plane.  We  then  count' 
cells  djxjdy  within  the  "reflection  twice.  Other  cells  in  the  upper  half  plane 
count  once  and  represent  electrons  which  reach  the  anode.  Still  other  cells 
lie  outside  W,  and  are  empty.  These  cells  lie  in  a kind  of  shadow  cast  by 
equations  (2.  3)  on  the  wings  of  the  (xy)-plane. 

The  virtual  cathode  lies,  of  course,  at  the  radius  r^  where  dijj/ds  is 
zero.  The  three  sections  in  Fig.  2-1  represent  W for  three  choices  of 
radius  r(or  s):  s >s>-l,8  = 8 and  s^*- s We  now  undertake  a more 

detailed  discussion  of  the  first  sketch  in  order  to  evaluate  S(r). 

In  Fig.  2- la,  W is  bounded  near  the  origin  of  the  (x,y)-plane  by  the 
y-axis,  and  in  the  wings  beyond  y^  by  dashed  curves,  segments  of  hyperbolae. 
The  hyperbolae  intersect  the  y-axis  at  y^  and  have  the  equations: 

X = 0 , X ^ = (s-l)y^  ; y^  = (2.  13) 


The  dashed  elliptic  segment  m represents  particles  which  just  reach  the 
virtual  cathode,  and  the  region  under  it  includes  most  of  the  cells  which 
count  twice  in  the  integral  S(ir).  The  segment  intersects  the  y-axis  at 
y^  and  has  the  equations: 


x^  + (1 


(2.  14) 


The  dashed  curves  and  equations  (2.  13)  and  (2.  14)  represent  roughly 
the  two  observations  made  above  on  the  boundaries  of  W.  W'e  note  that 
the  distribution  function  exp  ■|-(x^+y^+ij;)|  is  bell-shaped  and  centered  at 
the  origin  of  the  (x,y)-plane.  It  falls  off  rapidly  away  from  the  origin,  so 
the  part  of  W under  the  curve  m in  Fig.  2- la  is  the  most  important.  The 
excluded  wings  of  the  (xy)-plane  represent  very  small  values  of 
exp  |-(x“+y")J  , and  S(r)  retains  a value  near  exp  (-»j>)  as  long  as  j y^j  is 
larger  than  one. 


r 
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The  true  boundaries  of  W and  of  the  region  of  double  cells  are  not 

the  dashKu  curves  between  points  v-  and  y, , y - and  x , but  the  solid  lines 

* j 14  s 

running  diagonally  up  from  y^.  To  explain  the  solid  lines  we  introduce 
the  notion  of  an  operator  P(s)  which  can  transform  the  plane  of  emission 
velocities  into  the  plauie  (x,y)  of  velocities  at  the  radius  s:  (x,y)  = P(s)  (x,y). 
Under  action  of  P(s),  we  have  seen,  the  y-axis  on  the  (x,y)-plane  maps  into 
the  hyperbolic  segments  (2.  13)  ending  on  the  (x,y)-plane  at  y^^.  P(s)  is  thus 
an  abbreviation  for  equations  (2.3). 

We  consider  also  an  inverse  mapping  P ^(s*)(x,y)  = (x,y),  which 
carries  features  of  an  (x,y)-plane,  at  a radius  that  v/e  distineuish  from  s 
witii  a star,  back  on  the  (x,y)-plane.  In  particular  the  y-axis  maps  under 
P ^(s  into  an  elliptic  segment  with  the  equalionc  { ;};  O)'- 

x = 0,  + y^(l  - -^1  = iji(6+)  ; ■ (2-15) 

For  a series  of  values  of  s*  less  than  s , P”^(s*)  represents  the 

condition  x(a’*)  = 0 by  a family  of  elliptic  segments  with  a smooth  envelope 

that  extends  on  the  (x,y)-plane  from  y^  = + Vd  ijj Ms  , x = 0,  to  y = 0, 

“ V The  cells  dy  dx  lying  outside  this  envelope  represent  electrons 

which  will  cross  r and  reach  the  anode.  The  electrons  whose  emission 
m 

velocities  lie  inside  this  envelope  will  return  to  the  cathode.  The  direct 
mapping  P(s)  of  this  envelope  on  the  (xy)-plane  in  Fig.  2- la  yields  the 
solid  lines  through  y^.  The  points  y^  represent  points  y^  where  the  par- 
ticular ellipse  for  s = s*  touches  the  envelope.  The  points  y^  are  equal 
to  y^/Vs  . 

Since  the  mapping  P is  a continuous  function  of  s,  each  point  (x,y) 
i.r«u.e5  a curve  through  a succession  of  (x,y)-planes.  The  vector  ^tangent 
to  such  a curve  can  be  defined  almost  everywhere  on  the  (x,y) -plane  with 
equations  (2.  3).  Introducing  unit  vectors  x and  y along  the  x-  and  y-axes, 
and  partial  derivatives  to  imply  that  x and  y are  held  fixed,  we  write  t: 


t = 2 ^ + V ^ = X r4-»  - J_  -OX 
- = *8s^^»8  -^lix  2it  ds^  y 2j 


Then  the  intersection  y-  is  defined  thus: 


(2.  16) 


FiG.2-lFL0W  PATTERNS  AND  BCXiNDARY  OF  W 
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8£  = 0 . 7,  = +/^  . 8 <8 

08  -2  — Vd* 


(2.17) 


We  may  now  consider  briefly  the  other  sections  of  Fig.  2-1.  In  each 
one  the  vectors  t,  indicated  by  the  arrows  in  direction,  but  not  in  magnitude, 
give  the  pattern  of  flow  for  a particular  radius  s.  The  pattern  is  symmetric 
about  the  origin,  but  at  s = s^  all  along  the  x-axis  there  is  a change  from 
downward  to  upward  flow.  The  points  y^^  and  y^  have  the  same  significance 
in  all  three  sketches.  In  the  second  sketch  the  points  y^  have  closed  to- 
gether on  the  origin  and  the  region  of  double  cells  imder  mhas  vanished. 

The  points  also  tend  to  close  on  the  origin,  and  actually  meet  there  when 
»!>  is  zero  and  Hence  the  scale  of  the  third  sketch  represents  a 

considerable  enlargement  of  the  scale  used  in  the  first  two.  In  this  last 
sketch  the  region  above  the  origin  out  to  is  also  excluded  from  W c.nd  lies 


in 


O 


hadow  of  the  virtual  cathode. 


The  points  y^,  it  is  clear  from  their  definition,  never  close  together. 

A rough  approximation  to  W would  then  be  a rectangle  bounded  by  a hori- 
zontal line  through  x^  and  vertical  lines  through  y2-  A still  cruder  approxi- 
mation is  the  common  one  [2]  which  takes  for  W the  entire  half  plane  boimded 
below  by  X = x_.  This  yields  for  S the  following  formula,  in  which  the  plus 
sign  applies  to  the  interval  r < r<C*' 


m 


S(x)  = exp  [-v|i(r)l  -[1+erflxJJ  , . (2.18) 

Here  erf  is  the  error  function  or  probability  integral,  given,  for  instance, 
in  Peirce's  tables  [4]  . For  a planar  ^.athode  W has  no  wings  and  mis  in- 
deed the  line  through  parallel  tc  the  y-axis,  so  (2.  18)  is  an  exact  rela- 
tion, with  i|i  a fimction  of  the  cartesian  coordinate  perpendicular  to  the 
cathode.  To  find  S more  precisely  for  a given  cylindrical  cathode,  one 
may  use  numerical  methods,  or  semi-analytic  forms  for  an  imder-  and  an 
over-estimate  of  S,  as  Wheatcroft  does. 

The  net  outward  current  through  a cylindrical  surface,  concentric 
with  the  cathode  <ind  of  radius  r,  must  be  independent  of  r and  will  be  given 
by  an  integral  = jj  dl.  Here  dl  is  defined  in  (2.  9)  and  W.  is  the  region 

^i 
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V/  of  the  <x,y)-plane,  minuH  the  section,  if  any,  which  counts  twice  in  S. 
W.  and  di  must  change  with  radius  so  as  to  keep  constant.  A rough 
approximation  for  can  be  deduced  from  Fig.  2- lb: 


la  =?=  Iq{V«  er£|y-J  • exp  (-141)  + 1 - erf  |Vs  yjJ 


(2.19) 


This  gives  too  large  a value  for  I^,  because  it  includes  in  the  two 
triangular  regions  with  corners  at  y^,  yj^,  and  y^  = 0- 

If  the  anode  voltage  is  made  sufficiently  negative,  the  barrier  potential 

ijj  moves  radially  out  to  the  anode  end  disappears.  Then  if  and  s^  are 

substituted  in  (2.  14)  for  and  s^.  Fig.  2-la  will  represent  W for  all 

values  of  r,  r < r <r  . Schottky  studied  thermionic  curreits  flowing  against 
o ^ 

such  a retarding  potential,  and  (2.  19),  evaluated  for  s = s^,  is  his  formula 
for  I [1]  . It  still  includes  extra  triangular  regions  (y3,y^,y2)* 

Fig.  2-la  these  are  far  from  the  origin  and  introduce  small  error  into  1^^. 
Indeed,  if  jy^  | is  larger  than  2,  (2.  19)  reduces,  with  an  error  under  1 per 


cent,  to: 


Vs  exp  (-4») 


(2.20) 


Schottky' s formula,  we  note  parenthetically,  is  often  quoted  for  very  large 
s (or  very  small  cathode  radius)  [5,  page  80]  . tp.  this  case.  Vs  erf  jy^j 
approaches  2 V^i»  and  erfjVs  7i|  approaches  erf  but  these  form" 

do  not  apply  in  our  diode. 

In  the  next  chapter  we  give  results  of  a numerical  calculation  based 
on  the  true  boundaries  of  W,  the  solid  lines  in  Fig.  2-1.  To  complete  here 
the  formulation  which  will  be  used  in  this  calculation,  we  require  the 
potential  ^ to  satisfy  Poisson's  equation  (1.  16).  This  we  now  write  with 
one  initial  condition  and  normalizing  definitions  like  (2.  2)  and  (2.  11): 

[®*P  (o+2?)]  - S(^)  , i|>’  s ^ ^ ® 1 (2-  21) 


= 4»(^=0)  = 0 , exp(a)  = 


2el  r 
o c 

vh«  kT 
o 


I - 

> 4-‘  = lET  • 


(2.  22) 


I 

I 


■v/i'rwaf; 


TR201  -19- 

c . Integration  Methods  and  . 

In  this  section  we  use  the  formulation  developed  above  to  interpret 
some  data  obtained  from  our  diode.  With  Boltzmann's  approximation, 

S(^)  = exp  [->Jj(§)]  . one  can  integrate  (2.21)  and  obtain  the  results  in  Fig. 

2.2.  This  approximation  is  a good  one,  as  we  noted  following  equation 
(2.  14),  in  most  of  the  interval  r^<;;  ^rn‘  ^*'t^t>ducing  it,  we  put  (2.  21) 
in  the  form: 

^(f)^  = exp  (f)  , f = « + 2^  - 4>  , f = ^ (2.23) 

One  may  show  by  substitution  that  a solution  of  (2.21)  is: 
exp(f)  = b^  cos~^  (a  + ^§)  , or 

(5)  = 2B  + 2 in  cos  (a  + 1^)  - 2 In  cos  a (2.24) 

The  parameters  a and  b in  (2.  24)  depend  on  exp(a)  ana  on  f^,  the  sub- 
script c teferring  to  the  cathode  at  B = 0: 

b = cos  a [exp  (o/2)]  , = sin  a[exp  (ct/2)]  . (2.25) 

Here  exp  (a)  represents  the  parameter  of  emission  current  (2.  22),  and 
f^  represents  the  electric  field  at  the  cathode:  f'^  = 2 + (erE/kT)^.  We 
note  in  (2.24)  that  4i(B)  satisfies  - 0,  as  It  shovild.  The  second  condi- 
tion on  is  that  it  have  a maximum  at  B , or  that  f = 2.  These  bvo 
conditions  permit  one  to  determine  b as  an  eigenvalue  of  (2.  23),  azid  by 
further  algebra  one  obtains  the  results  in  Fig.  2-2. 

The  ordinate  in  Fig,  2-2  is  the  height  of  the  potential  barrier  4<^. 
and  the  abscissa  k describes  its  radial  position  r :k  =(s  -l)*(s  +1)  \ 

The  solid  lines  are  contours  of  constant  a.  If  they  had  been  drawn  in,  con- 
tours of  would  be  almost  parallel  to  the  contours  of  a;  actually  4^^  is 
represented  by  occasional  underlined  numbers  on  the  a-contours.  The  light 
dashed  contours  I^^A^  were  calculated  by  (2.  20),  evaluated  at  s^.  Hence  they 
have  the  triangulj.r  inaccuracy  mentioned  for  (2.  19).  However,  setting 
erf  = 1.  as  we  did  in  (2. 19)  to  obtain  (2.  20),  introduces  additional  errors 
in  I^Aq  only  in  the  region  of  heavy  dashed  contours  A = erf  ly . 
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Since  S(g)  must  lie  between  exp  (-i|i)  and  exp  (-Jj)  in  the  interval 
upperbound  in  a to  the  error  introduced  by  the  Boltzmann 
approximation  is  In  2 = 0.  7.  Thic  is  sbc-ul  half  the  a-interval  used  in 
Fig.  2-2.  The  position  of  each  a-contour  relative  to  its  neighbors  is 
much  more  accurate.  The  points  numbered  1 through  5 near  the  origin, 

^m  ~ ‘kn  ~ ^ calculation  for  a = 10  based  on  (2.  18),  and  the 

position  of  the  dashed  line  joining  them,  relative  to  the  solid  contour  a = 

10,  indicates  the  accuracy  of  this  contour.  Computations  to  improve  Fig. 
2-2,  although  it  applies  to  any  cylindrical  diode  and  has  some  general 
interest,  do  not  seem  justified  while  a larger  uncertainty  remains  in  the 

I 

experimental  values  of  a(T). 

Figure  2-3  presents,  in  the  line  through  the  points  x,  our  data  on  the 
function  a(T).  The  points  are  based  on  extrapolated  from  curves  of 
temperature  limited  current,  and  on  T from  curves  of  current  in  a retard- 
ing field.  Figures  5-3  and  5-4  of  T.P..  185.  The  parameter  in  Fig.  2-3 
is  y^,  the  voltage  applied  to  the  cathode  heater.  The  extrapolations  to  1^ 
above  940°K(V^  =12)  are  uncertain;  and  values  of  T below  820 °(V,  = 8.  5) 
are  merely  estimates,  so  the  points  x become  lines  indicating  a bracket 

on  a(T).  For  instance,  the  point  V,  = 12v,  with  a = 9.  55  and  T = 94C°K 

2 * ^ 

CO 'responds  to  0.89  ma  per  cm  , emission  current  density  and  0.63  watts 

per  cm  , heater  power.  In  terms  of  these  quantities  one  can  find  the 

point  on  a conventional  power-emission  chart  ^6,  Fig.  11,  page  83]  . The 

val'ues  appear  reasonable. 

> 

Now  is  measured  with  high  anode  voltages  V^,  and  T with  incre- 
ments AV  , so  the  cathode  and  anode  work  functions , V and  V , scarcely 
a c w •' 

alfect  I^,  T or  the  points  a(T).  (One  may  perhaps  question  the  combining 
of  high  and  low  voltage  pathode  characteristics  in  a single  parameter.  ) 

The  line  Y in  Fig.  2-3,  on  the  other  hand,  and  the  curves  in  Fig.  2-4  use 
only  data  on  the  diode  current  in  a retarding  field  and  yield  estimates  of 
V and  V -V  . For  Y,  data  from  the  linear  portion  of  the  curves  in  Fig. 
5-4,  T.R.  185,  were  used  (chiefly  the  intersection  V^(T)  of  each  curve 
■with  the  line  1=1  microampere).  The  formula  for  Y comes  by  eliminating 

3i 

between  Schottky's  equation  (2.  20)  and  Richardson's  equation  [7,  pages 
25  to  47]  : 
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i 


- eV 
c 

If  the  coefficient  A here  is  independent  of  T*  then  Y is  a.  linear  function  of 


1 = In  yx)  + 2 In  1 X ) + |x  V^l  , (2.  27) 

= -x(V^-Vy  + In  [AV'i^  (e/k)^]  , 

- a -I  |x  j + In  |x  I - 10.  96  ; x = (e/kT) 

Here  the  number  10.96  represents 


The  points  Y(T,V^)  fall  on  portions  of  two  lines  which  meet  in  = 

24  volts.  One  finds  A in  amps/cm^/°K  from  the  intersection  of  Y(x)  (ex- 
tended) with  X = 0,  and  V -V  in  volts  from  the  slope  dY/dx.  From  the 

- 2 ^ -3 

upper  line,  we  find  10  >A>10  and  1.5>  V -V  >1.3,  both  reason- 

c w ^ _ 

able  values;  but  from  the  lower  line  (V£  = 24  to  10  volts)  A = 27  x 10  and 

V - V =3.0,  quite  impossible  values  (5,  page  109]  . The  explanation 
c 

may  lie  in  heat  losses  through  the  ends  of  the  cathode  sleeve.  These 
losses  tend  to  overtake  radiation  losses  at  lew  lvalues  of  T and  may  give 
the  cathode  a non-uniform  temperature  distribution.  From  comparing 
o and  Y at  = 12,  we  find  in  (2.  27),  = 0,  26  volts. 

Figure  2-'4b  reproduces  portions  of  four  curves  of  diode  current 
against  a retarding  voltage  (T.R.  185,  Fig.  5-4).  These  curves  have 
been  tailored  in  two  ways:  The  slope  of  each  curve  in  its  linear  region 
was  reduced  to  unity  by  choice  of  T in  ijj*  = e(V^+V^y  kT  , then  each 
curve  was  made  asymptotic  to  one  line,  the  Schottky  line  (2.  20),  by  choice 
of  appropriate  1^.  The  selected  in  this  manner  proves  reasonably  con- 
sistent with  data  on  a in  Fig.  2-3.  but  a correct  selection  depends  criti- 

callv  on  a correct  value  for  V 
■'  w 

Normalization  of  and  with  respect  to  T^  and  has  not  changed 
the  shape  of  the  curves  in  Fig.  2-4b,  however.  Their  topology  resembles 


= AT^  exp  I 


' kT 


(2. 26) 
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that  of  the  curves  in  2-4a,  and  these  curves  were  computed  Irom  a-contours 
in  Fig.  2-2.  Fach  curve  in  Fig.  2-4a  leaves  the  common  asymptote,  the 
Schottky  line,  when  the  virtual  cathode  lies  just  in  front  of  the  anode.  As 
the  retarding  voltage  i|<  is  reduced  further,  each  curve  in  Fig.  2-4a  may 
be  considered  the  locus  of  a point  U ( kp  ),  and  a corresponding  point  Q ’’‘(di.) 
moves  along  a corresponding  a-contour  in  Fig.  2-2  from  k = 0.675  toward 

a 

the  cathode,  k = 0.  As  Q*(ilj  ) crosses  contours  of  1 fl  in  the  direction  of 

a a o 

increasing  1 , so  Q (4<  ) moves  up  a curve  in  Fig.  2-4a  toward  higher  values 
of  In  b'lt  ^a‘  '^a^  rises  more  slowly  now  than  it  did  on  the  Schottky  line . 

(This  mental  experiment  with  provides  a key  to  the  meaning  of  Fig.  2-2.  ) 
In  Fig.  2-4a  the  curves  end  at  = 0 and  the  dot-dash  line  in  Fig.  2-2  shows 
how  far  Q*  has  come  down  the  a-contour.  Thus  the  line  i|*  = 0 in  Fig.  2-2 

d d 

represents  the  position  of  the  virtual  cathode,  when  anode  and  cathode  are  at 
the  same  potential. 

The  computations  for  Fig.  ?.-4a  involved  matching  the  solution  of 

Poisson's  equation  (2.  24)  for  r^  to  the  Langmuir  solution  (1. 23) 

for  the  interval  r <1  r^  r . We  conclude  this  section  with  a brief  examina- 
m — a 

tion  of  Langmuir's  series  as  a solution  of  (2.  21)  for  r>  r^.  Beyond  the 
virtual  cathode,  W takes  the  shape  shown  in  Fig.  2-lc,  and  S takes  the  form 


oo 

S(s)  = — ^ f dy  [1  - erf  x(y)  ] exp  {-(y^  + 4<)]  (2.28) 

-00 


Here  x(y)  represents  the  solid  curve  between  points  y^  and  x^  in  Fig.  2-lc, 
plus  dashed  hyperbolic  segments  out  toward  large  x^,so  clearly  x(y)>x^. 

If  x^  is  larger  than  2,  we  may  replace  the  error  function  in  (2.  28)  by  its 
asymptotic  form: 


= 2V 


exp  + x^  (y)]| 

x(7) 


(2.29) 


if  v/e  replace  the  solid  curve  bordering  W by  the  dashed  hyperbolic  segment 
m,  x(y)  along  it  is  given  by  (2.  14)  and: 
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exp 


I 

-oo 


exp 


[-  3iLt^}  dy 


,.—2,8 

•ji  + y (— 

m 


(2. 30) 


1) 


Because  f“y^l  -ppcars  in  art  exponent  in  (2.  ;^0)  only  small  values  of  jy  j 


—2  ® 

count  and  we  may  neglect  Y ^^e  denominator.  Performing  the 

integration  we  find; 


S(.)<  i .xp  (2.  31) 

Substituting  (2.31)  in  (2.21)  with  redefinitions  we  obtain: 

_ _ -1  _ _ _ 

4/'  exp  (a  + 5)  (4  V5T^  ) . 4)'=d4,/d5  ; (2.32) 

® • ‘i'm  ’ 


The  series  f,  (^)  in  (1.  23)  is  I."a.ngmuir's  solution  of  (2.  32),  with 

^ — .4/3 

Initial  conditions  atr  : »b=0,ib'=0>  and  di  - df^ 

m ^ L 


d””  ~ = 9 exp  g [ l6  Vir  , or  d = 0.  46 5 exp  (2a/3)  . 


V<i.  i 


34) 


The  equation  which  Langmuir  solved  was  originally  derived  for  a space 

charge  limited  cylindrical  diode  in  which  the  velocity  distribution  among 

electrons  was  ignored.  It  is,  we  find,  the  asymptotic  form  of  Poisson's 

cqt^ation,  including  a Maxwellian  distribution  of  emission  velocities,  wiifen 

4<  is  larger  than  four.  Licingmuir's  solution  is  not  necessarily  asymptotic, 

for  given  to  solutions  arising  from  a more  careful  consideration  of 

S(t)  near  r^.  It  was  used  to  represent  <K^)  in  the  interval 

with  the  condition  o = oc  + t + In  (1  A ) to  match  it  to  the  .Boltzmann 

a o 

solution  in  the  interval  0<C ^ ^ yielded  for  Fig.  Z-4».  lt« 
use  is  partially  justified  by  the  resemblances  between  the  curves  in 
Figs.  2-4a  and  2-4b, 


d.  Summary 

The  space  charge  inside  the  virtual  cathode  in  a thermionic  diode 
depends  on  the  emistion  current  and  temperature,  and  T,  and  on  one 


i 


1 

i 

! 
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i 

other  parameter,  such  as  or  The  dependence  of  this  third 

parameter  on  the  anode  voltage  can  be  found  if  the  whole  charge  dis- 
tribution even  to  the  anode  is  known.  depends  on  the  cathode  material, 
surface  area  and  temperature.  The  use  of  T as  a parameter  represents 
the  well-founded  assumction  that  the  emission  velocities  have  a Maxwellian 

( 

’ distribution. 

. if  the  particles  in  the  space  charge  do  not  interact  individually,  their 

entire  orbits  obey  equations  (2.3).  These  equations  cut  away  the  edges  of 
the  velocity  distribution  as  the  particles  move  out  from  the  cathode.  The 
shape  of  the  edges  determines  the  values  of  the  integrals  I(r)  and  S(r),  the 
^ radial  current  and  the  charge  density.  S(r)  determines  the  potential  function 

4i(r)  through  Poisson's  equation,  and  4^r)  through  (2.3)  in  turn  helps  determine 
■ S(r).  Space-chargeproblemsusually include  this  circle  of  self-consistency. 

Lev/  anode  voltages  or  negative  oaes  which  cause  the  virtual  cathode 

\ 

tn  move  out  toward  the  anode,  bring  the  whole  tube  in  a sense  within  the  cathode 
region.  Hence  they  are  used  to  examine  the  effects  of  cathode  parameters,  al- 
though may  also  be  found  with  high  anode  voltages,  k|j^  = " ‘I'a* 

high,  the  velocity  the  electrons  pick  up  in  the  radial  field  far  exceeds  the  dis- 
tribution of  thermal  velocities.  Then  Langmuir's  series  solution,  in  which  T 
and  are  no  longer  important  parameters,  describes  k)j(r). 
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Effects  of  Emission  Temperature  in  a Magnetic  Field 

a.  Introduction 

The  problem  of  describing  space  charge  in  front  of  an  -mitting  sur- 
face parallel  to  a steadv  magnetic  field  has  apparently  received  little  at- 
tention [1,  page  211]  . It  is  common  to  assume,  in  discussing  magnetron 
performance,  that  the  magnetic  field  has  little  or  no  effect  on  the  virtual 
cathode.  After  considering  the  effects  of  the  tangential  component  of  elec- 
tronic emission  velocities,  Twiss  in  1953,  it  is  true,  proposed  a drastic 
revision  of  conventional  magnetron  models,  the  so-called  single  and  double 
stream  models  for  the  whole  charge  distribution  from  cathode  to  anode 
[2,  pages  247  to  289]  . But  Twiss  too  failed  to  examine  details  of  the  charge 
distribution  close  to  the  cathode,  and  we  shall  not  review  his  results  here. 

The  magnetic  field  does  in  fact  modify  many  features  of  the  cathode 
region,  as  we  shall  find  on  extending  the  analysis  of  Chapter  U.  But  the 
largest  effects  of  the  magnetic  field  occur  in  the  Hull  cutoff  transition.  We 
present  a rough  calculation  of  the  thermal  effects  to  be  expected  in  this 
transition,  and  conclude  this  chapter  with  results  of  a numerical  calcula- 
tion, The  results  allow  us  to  contrast  some  features  of  the  cathode  region 
before  and  after  i.  ''poaing  a steady  magnetic  field. 

b.  The  Virtual  Cathode  in  a Magnetic  Field. 

To  extend  our  previous  analysis  we  adopt  the  same  definitions  s,  i|i, 

3c  and  ~y  (2.2),  representing  the  radial  coordinate,  the  potential  and  the 
normalized  electronic  velocities  in  a cylindrical  diode.  We  continue  to 
ignore  the  effects  of  collisions  and  to  suppress  the  z-coordinate . More- 
over, we  neglect  possible  effects  of  a magnetic  field  on  the  mechanism  of 
emission  and  on  properties  of  the  cathode  materials. 

The  steady  magnetic  field  3,  uniform  and  parallel  to  the  system  axis, 
z,  we  normalize  with  reepcci  to  kT: 


(3.  1) 
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{Ve  will  avoid  using  in  this  chapter  the  symbol  h as  the  axial  length  of  the 
cathode.)  This  magnetic  field  does  not  enter  i.he  energy  equation  (2.  I), but 
the  new  momentum  relation  is  {1. 15),  r (0  + u/2)  = constant.  Instead  of 
(2,  3)  WE  have: 

+ y^  + = X + y~’  , Vs  y + 8 h = y + h (3.2) 

On  eliminating  y between  these  relations,  we  find: 

x^  - x^  + X(»)  "•  (s-l)(y  + \^  h)^  (3.3) 

2 

Here^(s)  might  be  called  an  effective  potential, (s)  = + h (s-1),  but 

convenience  rather  than  any  physical  concept  justifies  its  use, 

0 

The  definitions  of  incremental  charge  density  dS  and  current  dl 
(2.4,  2.5),  involving  the  energy  relation,  remain  unchanged.  The  Jacobian 
of  a transformation  from  emission  velocity  variables  (x,y)  to  other  velocity 
variables  {x,y)  remains  ^ (2.t).  Hence  dS  and  dl  transform  as  before 
(2.9,  2.  12).  For  the  integral  S(r),  the  region  of  integration  W on  the  (x,y)- 
plane  depends  chiefly  on  one  observation:  for  any  electron  in  radial  out- 
ward motion  that  contributes  to  S(r),  its  velocity  x(r)  must  have  been 
greater  than  zero  from  the  cathode  continuously  to  the  radius  in  question. 

In  other  words,  except  perhaps  for  singular  cases,  the  orbit  equations  al- 
low each  electron  only  one  outer  turning  point.  This  observation  leads  us 
to  define  wings  of  the  (xy) -plane,  excluded  from  W as  before.  We  complete 
later  the  integration  over  electrons  falling  back  to  the  cathode. 

The  features  of  W sketched  in  Fig.  3-1  should  be  compared  withthoae 
sketched  in  Fig.  2-1.  The  hyperbolic  segments  (y^.y^**  • solid 

lines  (731^2^  boundary  of  W have  the  same  significance.  The  points 

themselves  are  given  by  formulae  like  (2.  13,  2.  17): 

= ; Y^=s{dy/d3)  ; sY3^  = (d)j/ds)j^^ 

X = 0 = (8x/^s)  = 0 (end  of  envelope)  (3.4) 

These  points  Y sY  + V®h  lie  symmetrically  about  y = -^h  rather  than 
about  the  originY  = 0,  and  has  replaced  i|j.  These  differences,  due  to 
h,  are  cieiiT  already  in  (3,  3).  Revereing  the  direction  of  B merely  changes 


b) 


0 


<0 


FIG.  3-1  FLOW  PATTERN  IN  W,  h>0 
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the  sign  o£  h,  and  '.auses  the  sketches  in  Fig.  3-1  to  be  reflected  about  the 
line  y ~ - The  flow  vector  t;  moreover,  is  now  vertical  above  y = -2 V^h 

rather  than  above  the  origin  (2.  16). 


il!  - 


l=ij—  -dJ  • 


(3.5) 


The  contour  bounding  W near  the  origin  of  the  nE,y) -plane  is  the  most 

important  one,  since  the  integrand  in  S has  the  bell  shape  mentioned  before, 

centered  above  the  origin.  At  the  virtual  cathode  where  dij/ds  vanishes  and 

a equaH  a ; ♦■he  boundary  touches  the  origin  from  one  side  only  (Fig.  3- lb). 

It  cuts  from  W more  of  the  (x,Y) -plane  at  the  origin  as  s increases  beyond 

s . However,  in  part  of  W the  downward  flow  continues  out  beyond  s to  a 

second  radius  s^  where  the  points  y^  meet  and  dyjds  vanishes  (Fig.  3- Id), 

Bevond  s all  downward  flow  ceases, but  between  s and  s some  electrons 
ri 

which  have  passed  the  virtual  cathode  are  reflected  by  the  magnetic  field  and 
fail  back  to  the  cathode  surface.  Hence  in  the  px'csence  of  a magnetic  field 
8 replaces  s as  the  top  of  the  cathode  barrier.  This  is  in  spite  of  the 
fact  that  s^  lies  in  a region  of  accelerating  electric  field,  given  by  E(r^)  = 
mu^r^y4e. 

The  solid  curve  CYji.'yj)  near  the  origin  is  the  only  boundary  of  W we 
shall  have  to  find  in  calculating  S(r).  The  equations  for  this  line  were  based 
on  mapping  a series  of  lines  ^^(s*)  = 0 on  to  the  (x.y)-plane  with  P (s*)  and 
remapping  the  envelope  of  the  family  of  elliptic  segments  thus  obtained  back 
on  to  a particular  (x.y)-plane  with  P(c)  (Chapter  II).  In  these  equations  then 
sift  fixed  and  variable,  but  lies  mside  of  a. 

sYg(s*)^  = (8^  d^yai)  . iT  ^(s*)  = y^(8  + ) - X(s) 

_ _ 2 2. 

The  tangent  to  the  solid  curve  at  dx^ /dy^,  is  real  as  long  as  dy^  /ds 

is  negative  and  d-y/ds  positive: 


/=—  v2  r.  8 dr  /dr 

L ds  / 


(3.7) 


Clearly,  to  find  the  boundary  of  W at  s,  we  must  know  the  function  y(s  ♦ ) and 
its  derivatives  in  the  w'hole  interval  !•<  8"^=  <C  ». 


I 
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The  corresponding  formulae  for  a planar  magnetron  are  rather  less 
simple  and  symmetric.  We  record  them  as  the  basis  for  one  further  com- 
ment on  s^:  Placing  the  planar  cathode  {the  reference  for  potential, 
at  a height  abuve  the  (y,z)~ plane  in  order  io  normalize  with  reference  to 
X , rather  than  to  the  anode  height  d,  we  introduce  planar  definitions  like  (2.  2): 


s = (x/xj  , 4»  = i 

The  uniform  magnetic  field  B,  parallel  to  the  z-2Lxis,  appears  in: 


m 


— , ./  m 


(3.8) 


0)  = eBAn 


The  orbit  equations  for  an  electron  after  emission  are: 
u^  + + 4*(b)  = , V + hs  = V + h . 


(3.9) 


(3.10) 


After  eliminating  v,  v/e  obtain  a relation  like  (3.  3),  especially  if  y in  (3.  3) 
is  much  less  than  Vs  h: 


r 1*. 


u - u“  + X=  (8“1)  h [2v  + sh]  , "X  2 + h^B-l)  . (3.11) 

From  the  resemblance  of  (3.  II)  to  (3.  3),  j y”  j<[<C  VTh,  it  follows  thg*  in  the 

important  region  near  the  origin  on  and  the  (u,v) -planes,  the  curves 

' _ 2 

bounding  W are  almost  identical.  The  planar  formulae  for  v^,  v^,  v^,  u^ 
and  t,  moreover,  resemble  the  cylindrical  ones  if  one  assumes  i{j  ^^h^(s-l) 
and  expands  X (3.4,  3.5,  and  3.6). 

Away  from  the  origin,  however,  (3. 11)  leads  to  a parabola,  and  the 
envelope  of  mappings  P ^(s*)(u=0)  is  an  unbounded  curve  tangent  toa  family 
of  parabolae  on  the  (u.v)-plane.  Hence  the  vanishing  of  (d')^/ds)  does  not  de- 
fine a release  height  in  a planar  magnetron.  We  conclude  that  one  can  over- 
emphasize the  importance  of  s^  in  a cylindrical  diode.  The  eclipse  of  the 
origin  by  the  curve  as  we  saw,  begins  at  s and  relatively  few  elec- 

trons  fall  back  to  the  cathode  from  radii  between  s and  s . 

m n 

As  one  increases  the  magnetic  field  in  a cylindrical  diode,  the  z^elease 
radius  moves  out  and  the  Hull  radius  contracts,  so  eventually  these  radii 
may  merge.  Before  this  happens,  however,  our  present  analysis  fails  be- 
cause it  ignores  the  effects  of  collisions  and  the  stable  orbits  within  the 


L - 
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cathode  region  into  wlilch  electrons  may  be  ]3ushed  and  trapped.  To  integrate 
S(r)  over  the  electrons  falling  back,  to  the  cathode  as  well,  we  nnust  distinguish 
two  cases.  For  a diode  in  the  regime  before  cutoff  the  release  radius  helps 
fix  the  boundary  around  cells  d(x)dy  in  W which  count  twice.  After  cutoff, the 
anode  radius  determines  this  boundary,  and  virtuaUy  all  cells  lying  within  W 
count  twice . 


jj 

*'•» 

i 

,'■* 

A 

A 

5 


c.  The  Cutoff  Transition. 

Before  cutoff,  in  other  words  before  8^  reaches  2w  (Fig.  1-1),  some 

electrons  which  were  released  from  r^  with  unfavorable  tangential  velocities 

fail  to  reach  the  a.acde.  They  are  then  lost  to  the  anode  current  and  in 

completing  their  paths  back  to  tlio  cathode  they  contribute  c.  double  shielding 

effect.  For  both  reasons  they  cause  to  decrease,  yet  the  analysis  in  Chapter 

1 took  no  account  of  such  electrons.  As  the  magnetic  field  increases  through 

its  critical  value,  an  ever-increasing  proportion  of  the  radial  current  through 

r misses  the  anode.  This  transition  is,  of  course,  the  Hull  cutoff.  We  now 
n 

wish  to  extend  the  W integration  clear  to  the  anode  and  find  the  ratio  I. A as 

ci  U 

a function  of  p.  We  consider  that  cutoff  begins  where  this  ratio  starts  to  fall 
below  one,  at  p = p.. 

Figure  3-2  may  clarify  the  argument.  The  first  sketch  represents  the 
plane  of  emission  velocities  (x,y):  the  solid  curve  is  the  envelope  of  the 
family  of  elliptic  segments  x(^)  = 0,  and  the  ellipse  x(s^)  = 0 touches  it  at 
y = h.  The  second  sketch  represents  the  velocity  distribution  at  the  release 
radius  s^:  the  solid  curve  inside  ~y^  represents  the  envelope,  and  the  y-axis 
is  the  dashed  ellipse  in  the  previous  sketch.  The  third  sketchis  an  enlarge- 
ment of  the  area  near  the  origin  in  the  second  one.  The  fourth  sketch  is  of 
this  same  area,  referred  to  axes  (u,v)  which  are  aligned  with  the  boundary  of 
W near  the  origin,  and  hence  lie  at  some  angle  y to  the  (xy)-axes. 

In  each  sketch  the  dot-dash  curves  represent  the  conditionals  ) = 0. 

In  the  fir?t  and  second  ones  a single  such  curve  appears,  an  elliptic  seg- 
ment, with  the  equations: 

x^  + (1  - ^)(y  + h)  =X(»a^  " h^(e^-l-8p)  in  3-2a  , (.3.12) 

x^  + (l  - s /s  )(y  = 3-2b  . (3.13) 

n ok  IX  A II 
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This  curve  passes  through  W only  if  i»  larger  than 

definition  of  p.  is; 

(Apparently,  in  a planar  magnetron  like  s^  moves  to  infinity,  and  the 
smallest  magnetic  field  reflects  some  electrons  from  all  heights  x back  to 
the  cathode.)  When  p equals  p^,  we  ass\ime  that  s^  still  lies  tairly  close  to 
the  Cii'.thode.  The  second  derivatives  of  i]j  and  X must  vanish  between  s^  and 
8 , and  both  curves  must  have  inflection  points.  In  fact,  while  |«|>|  beyond 
8^  is  a monotomic  increasing  function  of  s,X(®)  behaves  more  like  a poly- 
nomial cubic  in  8. 

In  the  third  and  the  fourth  sketches  three  dot -dash  curves  represent 
a corner  of  the  elliptic  segment  (3.  13)  for  three  successive  values  of  p. 

Now  in  fact  changing  p will  change  the  whole  function  x(s)  and  with  it  the 
solid  line  bounding  W,  but  we  have  to  assume  these  changes  are  small  com- 
pared to  the  movement  of  the  curve  ^8^^)  = 0*  This  means  we  regard  1^ 
as  fixed  during  the  cutoff  transition  and  ignore  the  extra  shielding  of  the 
virtual  cathode  by  electrons  whose  velocities  at  place  them  inside  the 

dot-dash  curve  x(s^)  = C. 

«> 

From  Fig.  3-2c  the  current  increment  dl(u,v)  and  from  Fig.  3-2d 
the  integral  1^  over  the  half  plane  can  be  written  by  inspection; 

dl  = e*P  f-(4<j^+'i^4-v^)j  [cos  y d(u^)dv  + sin  ydud(v^)]  ; (3.15) 

I = I i\/b  exp  +u  ^)]  cos  w . 

n on  L n o J ■ 

Here  u^  is  the  distance  along  the  u-axis  from  the  origin  to  the  edge  of  W^. 

The  current  I has  the  same  increments  as  1 , but  the  area  of  inte- 
a n’ 

qration  W is  reduced  from  W by  removing  the  cells  under  the  dot-dash 

^ A 1- 

curve.  This  curve  cuts  the  v-axis  at  v . We  take  v as  the  v-coordinate 

n n 

also  of  the  intersection  of  the  dot -dash  curve  with  the  s;>lid  edge  of 

As  p decreases,  the  ellipse  x(s^)  = 0 grows  and  v^  sweeps  across  the  origin. 

Our  problem  is  to  make  explicit  the  dependence  of  I on  p from,  the  formula: 
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Iq  r 2 

~ exp  ( -4<n)  I )cos  v ( 1 + erf  j ) + 

exp  (-v/)sin  Y (1  -erfu^)]  +\^exp  (-4,^)(1  -erf|7J)]  (3.16) 

less  than  zero,  we  choose  the  plus  sign  before  erf  |v  j,  etc.)  The 

liFt  *rT*T*  in  (3.  16}  represents  tKe  in  Fif>  «i 

P(Sj^)(x,y)(Sn)  ~ (^•yH8a.^*  velocity  variables  at  the  anode,  maps  the  dot- 
dash  curve  into  the  y(8^)-axis  and  yields  this  term. 

The  formula  used  to  compute  the  two  curves  in  Fig.  3-3  contains  only 

exp  (-y^'O  I 


the  expreaBlons  K and  y^;  it  is: 


^a  1 

r~"  I 

n 


(1  + erf  |y^|)  + 


»V7(i-K)j 


(3.17) 


To  derive  it  from  (3. 15)  and  (3. 16)  we  first  assume,  in  view  of  our  ignorance 
ox  iae  exact  shape  oi  '*7^,  that  y and  u^  are  zero  and  v equals  Y • Then  we 
have: 

^a  1 f,  , . . — i/®^ 

(3. 18) 


a 1 


f iP 

(1  + erf  |y^|)  + exp  d - «rf  | yj)  | 


^ 2 

We  define  y^  by  setting  x = 0 in  (3.  13)  and  dropping  4>(s^)  as  in  (3.  14): 


y £^^~hf-  1 + 

XI 


a 


s -a 
a 


:1 


(3.  19) 


The  point  representing  y on  the  (■x,y){s^)  plane  is  y",  : 


y = yT~  h 
^a  a 


(3.20) 


Whenever  y^  is  near  zero,  y^  is  less  than  zero  and  | y^  |is  larger  than  one 
(K  = s /s  < 0.  2),  so  we  can  write: 


1 -erf|y  I = 


I - [exp  (-y  ^)]  (xy  ^) 
' I.  a j a 


(3.21) 


Finally,  from  the  energy  relation  (3.  2)  with  x(s  ) = 0 = x(s  ),  we  have- 

a n 
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+(3^^  (3.22) 

Thus  we  arrive  at  (3. 17),  giving  a*  a function  of  p and  h.  The  center 

of  the  cutoff  transition  occurs  whenV  is  z^ero,  and  this  defines  p <C  P. : 

'n  *^c  1 

Pc  = “Ti—  = ST—  ®a>>  ®n  = " • 

A 


In  general,  the  two  curves  in  Fig.  3-3  show  how  much  a distribution 
in  emission  velocities  would  tend  to  round  off  a transition  through  cutoff,  if 
the  current  through  the  cathode  barrier  could  be  fixed  while  the  magnetic 
field  increased  or  p decreased.  Actiially  1^  is  not  fixed.  One  might  expect 
1 to  fall  by  half  in  the  transition,  since  electrons  which  double  back  from 
anode  to  cathode  roughly  double  their  contribution  to  shielding  the  virtual 
cathode. 


The  curves  in  Fig.  3-3  should  be  mviltiplied  by  before  standing 

comparison  with  ihe  data  in  Fig.  3-4.  If  also  decreased  in  the  cutoff 

transition,  the  calculated  curves  would  be  sharper  than  3-3,  but  we  find 
the  data  in  Fig.  3-4  shows  rather  a less  abrupt  cutoff.  This  suggests  that 
another  transport  niechanism,  a55c..?i~tcd  perhaps  with  collisions  between 
charges  moving  radially  in  and  out,  affects  the  shape  of  the  transition. 


The  data  in  Fig.  3-4  for  the  regime  before  cutoff,  from  p = 10  to 
Pj^  = 0.55  is  in  qualitative  agreement  with  the  transit  time  effects  cal- 
culated for  Fig.  l-l.  On  the  other  side  of  p^,  the  curve  = 100  in  Fig. 
3"4  is  more  abrupt  than  the  curves  = 2.5;  indeed,  the  angle  between 
them  is  larger  than  between  the  calculated  curves  in  Fig.  3-3.  Of  the  two 
curves  V = 2.5,  the  one  taken  with  the  hotter  cathode  (T  = 1230°K,  V,  = 30) 

& Cl 

as  '.VO  would  expect  is  slightly  more  rounded  near  p = p.  than  the  other  (T^  = 

1130°,  y,  = 24  voiis';.  To  make  the  cux'ves  V = 100  and  V =2.5  cross 
1 a a 

near  p^,  0.75  volts  was  added  to  the  potential  applied  to  the  tube  in  cal- 
culating p.  V/e  saw  in  Chapter  II  that  - V„.had  this  sign  and  perhaps 
tvdce  this  magnitude. 


d.  A Numerical  Example 


In  this  section  we  show  by  numerical  work  the  dependence  of  the  cathode 


FIGURE  3-3  FIGURE  3-4 
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region  on  a,  h and  The  calculations  carry  us  only  to  '§  = 0.  06  or  0.  P2 

cm  from  the  physical  cathode  in  our  diode.  The  calculation  for  the  case 

with  magnetic  field  stoos  then  between  s and  a , yet  to  carrv  it  further 

m n 

without  a consideration  of  trapping  orbits  and  collisions  seemed  unrealis- 
tic. The  assumption  that  all  regions  of  W count  twice,  that  the  Hull  radius 
is  less  than  the  anode  radius,  much  simplifies  the  calculation  for  the  cathode 
in  a magnetic  field.  In  the  case  without  magnetic  field  one  must  distinguish 
single  and  double  regions  on  the  (x,y)-plane,  so  the  integration  starts  at 
4<m>  moves  radially  both  outward  ard  in  toward  the  emitting  surface. 

To  give  a few  details,  we  integrate  Poisson's  equation  (2.21)  with  the 
initial  conditions  si.'  =0,  and  the  set  of  six  values  of  ib'  shown  in  Fig.  3-5 
for  h - 10,  and  five  values  of  ili’  for  h = 0.  The  values  a = 10  and  h = 10  are 
appropriate  to  our  diode  in  a magnetic  field  after  cutoff.  Figure  3-5  shows 
solid  curves  of  i{i  and  ij*'  for  h = 10,  and  dashed  curves  for  h = 0. 

Each  solid  curve  was  started  by  Wheatcroft's  solution  (2.  24)  and  con- 
tinued in  this  form  until  S(^)  fell  appreciably  from  exp(-s|;),  the  value  as- 
sumed in  (2.  24) , The  transition  to  a numerical  method  at  actually 

occurs  when,  in  each  case,  y^  as  given  in  (3.4)  equals  1.5.  We  find  S is 
still  greater  than  0.95  exp(-v|;)  at^j.  From  "gj  we  carried  the  integration 
forward  in  steps  d^=  ^^i  “ ^i  i’  <3.  6)  at  each  step  to  plot  the  boundary 

of  W^.  To  evaluate  the  integral  of  the  distribution  exp  |^-(x^  + y^)J'  we  placed 
over  each  W.  a transparent  sheet  divided  in  0.2  x 0.2  squares,  each  labelled 
according  to  its  weight  in  the  distribution.  We  then  added  up  the  weights  of 
the  squares  within  W,  including  fractional  weights  for  squares  cut  by  the 
boundary.  The  total  was  exp(t|i.).  Conventional  integration  formulae  then 
gave  and  i|>^,  which  were  used  in  turn  to  compute  the  boundary  of  W.^j^. 

Figure  3-6  carries  further  the  comparison  of  the  solid  and  dashed 

curves.  We  see  there,  for  given  a and  tb'  , that  S occurs  earlier  and  vb 

is  lower  for  h = 10  than  for  h = 0.  Indeed,  J>  and  t , h = 0,  areaboutthe 

same  for  a given  initial  condition  di'  as  tb  and  t , h = 10,  for  ib'  +15. 

c ^m  -^m  ^c 

The  displacement,  15,  in  is  due  chiefly  to  our  assumption  thatwithh  = 10 
the  magnetron  diode  current  is  cutoff,  and  that  all  parts  of  W have  as  many 
falling  electrons  as  rising  ones.  Consequently  the  normalized  charge  density 
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S at  an  arbitrary  point  beyond  (actually  the  second  zero  of  »(<)  is  nearly 
twice  as  large  for  a given  when  h is  ten  a>>  it  is  when  h is  zero.  The 
dashed  lines  in  Fig.  3-b  repreeent 

In  integrating  the  dashed  curves  back  from  the  virtual  cathode  to^  = 0, 

we  had  to  estimate  the  pairs  of  initial  conditions  CH  4*  ) which  would  yield 

^ ’in  ' 

a = 10  at  ^ = 0.  Theae  estimates  resulted  in  slight  deviations  of  a from  10. 
The  deviations  range  from  C.02  for  curve  1 in  Fig.  3-5  to  0. 10  for  curve  5. 

e.  Review 

Figure  3-7a  shows  clearly  the  bell -like  shape  of  the  distribution  func- 
tion exp|^-(x  + y ) I . We  have  seen  how  the  equations  of  the  electronic 
orbit,  (2.  3)  or  (3.  2)>  cut  away  the  wings  of  this  distribution.  The  boundary 
of  the  cut  depends  on  the  shape  of  the  potential  function,  44^*)  or  'y^(s'*‘);  in- 
deed the  y'-coordinates  of  points  on  the  boundary  are  just  the  slopes  cf  cer- 
tain tangents  and  sectors  shown  in  Fig.  3-?b.  The  sketch  shows  clearly  the 

IL  2 ^ 2w 

progression  y^  ^ y^^  > {y with  h ^ 0,  the  same  progres- 
sion in  could  be  found  or.  a curve  Xi®* 

In  Chapters  II  and  111  we  have  concentrated  on  the  cathode  region.  In 
Chapters  IV  and  V we  take  a larger  view  and  regard  this  region  as  a sur- 
face layer  on  the  space  between  physical  cathode  and  anode  cylinders.  As 
the  source  of  electrons  for  the  interval  s<^  s^,  which  we  may  call  the 

interaction  region,  the  cathode  layer  determines  the  number  and  velocity 
distribution  of  the  electrons  that  enter  th<-  interaction  region.  In  this  aense 

a,  h,  and  d<'  are  independent  parameters,  but  in  another  sense  a > ® >®  > 

c m m n 

and  the  other  features  of  the  surface  layer  depend  on  the  anode  voltage  V^, 
and  on  the  distribution  of  space  charge  within  the  interaction  region. 

In  the  diode  before  cutoff,  and  even  rather  close  to  the  cutoff  transition, 
the  Langmuir  potential  function  (1.23,  2.32)  gives  4>(®*)  through  the  interac- 
tion region,  for  positive  and  large  enough.  This  function,  fitted  in  the 
catnode  layer  to  curve  2 of  Fig.  3-5a,  appears  as  in  Fig.  '^-7c.  The  curves 
S + h^(s-l)  in  this  figure  represent  two  cases:  for  h = 6.6,Xj__j  corre- 
sponds to  p p^  where  thermal  effects  on  the  cutoff  transition  begin,  and  for 
h = 7.5,  corresponds  to  cutoff,  p = p^.  The  radii  s^  and  are 
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FIG  3-6  COMPARISON  OF  CASES  h = 0,  h = !0,  FOR  SAME 
EMISSION,  a =!0,  AND  VARI0L6  ELECTRIC  FIELDS 
AT  THE  CATHODE 
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iitdicated  for  Xlj*  curve  has  a Sccoiid  maximum  at  s^,  a semicubic 
shape,  and  an  inflection  point  between  1 aud  ?■.  although  n.  the  charge  density, 

I 

does  not  vanish  in  this  interval  of  s.  One  may  also  compare  (1.13)  of  T.R. 

7. 

185,  the  expression  for  u , with  this  expression  for  the  inflection  point; 


er 

c 

4r^kT 


rZE 


*■  r £_■* 


(3.24) 


The  dashed  curves  in  Fig.  3-7c  represent  a quadratic  potential  func- 
tion 4<q  and  a corrcspcnding  Xq  curve.  At  cutoff,  the  correct  curve  for  X 
is  probably  closer  1°  than  to  because  of  collision  interactions. 

These,  we  will  argue  in  Chapter  V,  proviae  a charge  transport  mechanism 
which  seems  to  move  the  center  oi  gravity  of  the  whole  charge  cloud  out 
away  from  the  cathode  layer. 
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IV 

The  Magnetron  Diode  After  Cutoff 

a.  Introduction 

Most  models  for  the  distribution  of  space  charge  within  a magnetron 

diode  after  cutoff  provide  at  the  Hull  radius  s„  an  outer  boundary  for  the 

rl 

charge  cloud.  In  this  chapter  we  consider  three  well  known  models  of  the 
bounded  cloud.  After  a review  of  the  probe  signal  evidence,  we  reject  these 
models  in  favor  of  an  "unbounded”  cloud.  We  conclude  that  the  space  charge 
entirely  fills  the  interaction  space  in  the  diode  and  formulate  the  potential 
function  for  this  model. 

b.  Models  of  a Bounded  Space  Charge  Cloud 


Hull,  as  we  noted,  based  his  ideas  in  1921  [1]  on  Langmuir's  model  of 
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a cylindrical  diode  without  magnetic  field  since  formulae  in  T.R.  185 

occur  often  in  this  chapter,  we  indicate  references  to  that  report  by  an 
asterisk).  V/ith  the  energy  and  mcrr.cnturr.  relations  (1.4*,  2.  1),  Hull  ex- 
plained the  sharp  cutoff  of  anode  current  at  the  critical  ratio  p^,  of  anode 
voltage  V to  the  square  of  the  magnetic  field.  If  these  relations  apply  in 
the  same  way  after  cutoff,  if  virtually  all  electrons  within  the  space-charge 
cloud  retain  the  energies  and  the  canonical  momenta  that  they  have  onleaving 
the  cathode,  then  the  cloud  consists  of. electrons  in  emission  orbits.  Such  a 
cloud  will  be  bounded  by  a radius  s^^  equal  to  the  outer  turning  radinr  s^  of 
the  typical  emission  orbit. 

The  potential  at  = s^  (1. 7*  ) in  this  case  must  pass  smoothly  into 
a iogaritlimie  potential  appropriate  to  the  charge -free  cylindrical  region 
beyond  Sj^: 

2 


*(.) . - 


2 (^)  1 = — 5 


m CO  r. 


(8b-l)‘ 


(4.  1) 


'*>(8)  =<*>(Sj^)  + K(Sjj)  In  , Sjj<  s<  s^  ; K(s)  .5  . (4.2) 

H 


If  the  quadratic  expansion  for  potential  (1.  14*)  represents  the  potential 
throughout  the  radial  extent  of  the  cloud,  we  have: 

«>{8)  = {^)[B  + C(8+D]  . K{s)  = ^[C8  + 5i£]  . l<8<8j^  (4.3) 

If  the  electric  field  in  this  case  vanishes  at  the  cathode,  K(l)  = 0,  then  with 
(4.1)  we  find  the  constants  B and  C:  C = 1,  B = -2C.  The  resulting  potential, 
4>(s)  = (s-l)^/is,  was  proposed  by  Brillouin  in  1941  [2]  for  the  region  within 
the  cloud.  It  requires  the  canonical  angular  momentum  P to  equal  one 
throughout  the  cloud;  hence  the  circumferential  velocity  Y vanishes  at  the 
C^taaCde,  dLU  it  S^lOU  average  emission  component  is  zero: 

P = fQ  + '4l  = i ; Y = , Y(r  ) = 0 . (4.4) 

2 

Let  us  examine  briefly  Brillouin' s famous  model:  Setting  1 - C = 0 = u 

implies  the  absence  of  radial  motion.  We  note  from  (1.  15*)  that  the  state 
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u <:,  0 is  unstable,  for  if  n is  negative,  s becomes  a hyperbolic  function  of 
time.  Then  the  electrons  soon  hit  a surface  of  iLe  tube-  or  ot  leave 

the  region  where  u^  is  less  than  zero.  Since  u^  represents  in  part  the  charge 
density  (i.  13*),  a charge  cloud  with  negative  u^  soon  thins  out  until  u be- 
comes positive.  Brillouin's  model  lies  at  the  edge  of  this  unstable  condition, 
and  Brillouin's  cloud  has  the  maximum  density  permitted  in  a given  magnetic 
field.  Since  in  his  model  the  electrons  have  no  radial  motion,  they  move 
around  the  cathode  as  a single  slipping  stream.  He  devotes  much  of  his 
most  recent  paper  on  the  magnetron  [3j  to  showing  how,  as  one  slowly  raises 
V from  zero,  a single  cloud  in  such  motion  can  form. 

A 

Brillouin  there  also  states  and  proves  a theorem  that  the  double -stream 
model,  a rival  model  proposed  by  Page  and  Adams  in  1946  [4]  , is  only  pos- 
sible for  Vs jj  ^ 2 . 27  3 . However,  this  theorem  and  both  models  depend  on 
the  assumption  that  the  thermal  distribution  of  emission  velocities  may  be 
ignored.  Twiss  has  shown  [5]  that  this  assumption  and  the  uncritical  use 
of  Llewellyn's  equation  lead  to  singular  rather  than  general  results.  Nothing 
significant  in  magnetron  data  occurs  at  2.27  3. 

Let  us  also  examine  Ihe  model  sf  Page  and  Adams;  In  it  the  space 

charge  consists  of  two  equal  and  opposite  radial  currents,  and  ~I]^>  and 

also  flows  as  one  circumfereiitial  current  around  the  cathode.  The  radial 

motion  requires  a potential  to  which  the  quadratic  form  (4.  3)  cannot  come 

2 

very  close,  but  one  can  show  that  u lies  in  the  interval  0.  15  to  0.075.  Hence 

the  charge  density  in  this  model  is  not  very  different  from  Brillouin's.  The 

2 2 
probe  signal  data  includes  values  of  u from  one  to  0.  25,  but  no  values  of  u 

were  found  closer  to  zero.  Still,  one  must  add,  such  resonances  m.ight  be 

too  broad  to  ?hcv.’  at  all. 

Twiss  [5]  recently  described  a multiple  stream  model,  not  too  different 
from  the  model  of  Page  and  Adams.  In  his  set  of  initial  assumptions,  how- 
ever, he  included  a thermal  distribution  of  emissibn  velocities.  In  Chapter 
III  this  led  us  to  define  a release  radius  in  a region  of  accelerating  field. 
Disregarding  other  effects  of  emission  velocities,  we  can  easily  include  the 
electric  field  at  the  virtual  cathode,  K(l),  as  an  adjustable  parameter  in  the 
double  stream  model.  Ignoring  many  elegant  aspects  of  the  treatment  by 
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Twiss,  we  will  develop  this  model  in  the  rest  of  this  section.  The  reasons 
we  give  for  rejecting  it  will  apply  as  well  to  his  treatment. 

T:j  find  u in  the  modified  double  stream  model,  we  define  M (1.  IS) 
in  terms  of  I,,  the  outward  component  of  iTcttil&l  Ciix^cikty  AAid  tllCn  uCVclOp 
a form  like  (1.  19): 


M = el,  (Zirhe  m w^r  ^ 
1 o c 


(4.5) 


\ Vl^\  1 . r (1*^®  - Ins*)  ds* 

4>(s)  = K(l)  Ins  + M J 

1 Vss*  4»(a*)  - (b*-1)' 


, K S ■<  B, 


(4.6) 


We  may  develop  <^(s). as  a function  of  the  parameters  K,  M and  s^,  by 
iteration  over  the  charge-free  potential 4>^(s*  ),  and  express  M and  K in 
terms  of  s„  (or  s.  ) and  a new  parameter  m: 


4>(Sj^)  In  s* 


(l-m)*(Sjj) 


(4.7) 


_ M 


(Ins  - ln&»)  ds» 

'^1  Vso*  4*(s*)  - (s*-l)“ 


(4.8) 


Here<t(8  ) is  given  in  (4.  1),  and  we  extend ^(s)  in  the  interval  s,,<r.s-<6 
xi  XX  a 

with  (4.Z-).  Thus  we  find  p(m,s,  ) =^(s  ).  Knowing ^(s),  we  may  also 
2 D a 

evaluate  u (m,s^): 


= w T,  = w a 


(4.9) 


1 VSs'iKs)  - (s-l)* 


in  Fig.  4- la  we  plot  the  inverse  relations,  contours  of  m and  k on 
2 

the  (p,u  ) -plane.  The  quantity  m represents  the  field  at  the  cathode:  thus 
if  K(l)  = 0,  m = 1,  and  if  m = 0 then  K(  1),  for  a given  p,  has  its  maximum, 
a value  appropriate  to  the  cathode  field  in  a charge -free  coaxial  region. 
The  quantity  k represents  the  radial  extent  of  the  cloud  of  charge: 

®K  ■ ^ 

b c 


1 


(4.  10) 
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Although  we  made  only  one  iteration  over  for  the  *P(b)  used  in  the  (m,k)- 
grid.  Vv'e  helie-i/e  from  reault*  of  a,  direct  calculation  for  a planar  diode  that 
further  iterations  would  not  change  this  grid  much.  The  absorption  data 
of  Figr  4- 2+  interpreted  with  the  grid  on  Fig.  4- la,  imply  very  high  values 
of  K(l)  as  p approaches  zero.  The  data  do  not  suggest  the  line  K(l)  = = 0, 

m = 1,  appropriate  to  the  Page  and  Adams  model. 

The  parameter  m (or  M)  also  gives  the  charge  density  p in  the  cloud 

! T ^ 

(1.  18,  M = peb/4m€  Introducing  the  symbol  q for  the  average  (barred) 

space -charge  density  inside  Sj^,  we  write: 

s 


— Pe  f’  1 4Mu)  t 4irM 


(4.11) 


t. 


2 

From  this  one  obtains  the  q-grid  on  the  (p,u  )-plane.  Fig.  4-lb.  The  absorp- 
tion data  do  not  fall  along  any  contour  of  average  charge  density  normalized 
with  respect  to  However,  dividing  ^ by  p serves  to  renormalize  p,with 

respect  now  to  V^.  This  yields  the  contours  N = q/p  = pr^  /^a 
4-lc.  The  absorption  data  tend  to  follow  the  line  N = 1.5.  This  suggests 
that  the  anode  voltage  determines  p,  and  actually  led  us  to  try  the  quite  dif- 
ferent model  described  below. 

c.  Evidence  for  an  Unbounded  Charge  Cloud 

We  have  seen  that  modifying  the  double  stream  model  to  include  a 

finite  electric  field  at  the  cathode  is  not  enough  to  interpret  the  probe  signal 

data.  Let  us  then  consider  relaxing  the  energy  condition.  We  noted  that  the 

charge  cloud  is  bounded  by  s^^  in  any  model  in  which  individual  interactions 

do  not  disturb  the  electronic  orbits.  We  saw  that  s„  must  retreat  toward  the 

Jn 

cathode  as  p falls  from  p^  toward  zero.  However,  if  interactions  do  permit 
electrons  to  exchange  energy  and  momentum  at  a significant  rate,  then  the 
space-charga  cloud  may  extend  to  the  anode  for  all  values  of  p.  Hull  ap- 
parently had  such  a model  in  mind  in  19^4  when  in  au  abstract  [6]  he  sug- 
gested that  the  cloud  formed  "a  solid  shaft  of  rotating  electricity." 

There  is  little  experimental  evidence  for  a well-defiiled  outer  boundary 
to  the  charge  cloud.  Reverdin  [7]  four  years  ago  passed  a wide  beam  of 
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electrons  parallel  to  the  polar  axis  through  a grid,  the  magnetron  diode  and 
a second  grid.  On  a phosphor  screen  he  observed  the  beam,  the  shadow 
cast  by  these  grids,  and  uie  distortion  in  these  shadows  due  to  the  beam's 
having  passed  through  the  radial  field  in  the  magnetron  cloud.  He  found  at 
first  indications  of  a maximum  density  near  s , but  as  he  brought  his  cathode 
into  alignment  with  the  beam  and  the  magnetic  field,  the  density  became  a 
smooth  decreasing  function  of  radius.  Peterson  [8]  this  past  year  passed  a 
very  narrow  electron  beam  cixially  through  the  magnetron  cloud,  observed 
its  deflection,  and  concluded  there  was  considerable  charge  outBlde  the 
nominal  Hull  radius.  Knoll  and  Nedderman  [9]  have  used  excited  helium 
atoms  to  indicate  the  space -charge  distribution.  They  found  "as  the  anode 
voltage  is  lowered,  the  charge  density  decreases,  and  the  position  of  maxi- 
mum density  shifts  toward  the  anode." 

Several  experiments  have  detected  electron-  which  gained  great  energy 
within  the  interaction  space  and  were  not,  therefore,  in  emission  orbits.  For 
example,  we  recall  the  experiments  on  cathode  heating  by  back  bombardment 
[10]  . Wigdortschik  [11]  and  Linder  [12],  as  further  examples,  measured 
currents  to  "hats"  at  the  ends  of  the  diode  axis  when  these  hats  were  much 
more  negative  than  the  cathode,  hundreds  of  volts  move  negative  when  the 
anode  was  over  a thousand  volts  positive. 

Moreover,  the  electrons  collected  as  leakage  current  at  the  anode  are 

not  in  emission  orbits.  As  p falls,  the  anode  current  1 continues  to  flow; 

and  in  any  bounded  cloud  model,  we  must,  it  would  seem,  imagine  a mechanism 

that  knocks  some  electrons  out  of  the  cloud  and  clear  to  the  anode,  yet  does 

not  fill  the  stable  orbits  that  lie  like  traps  between  s,-  and  s . The  curves 

ft 

ir  the  four  square  sections  of  Fig.  4-2  show  1 over  a range  of  p after  cutoff. 

Each  section  gives  data  for  a fixed  V , and  within  each  section  the  cathode 

& 

temperature  is  a parameter.  The  coordinates,  p and  before 

in  Fig.  3-4,  here  are  plotted  on  logarithmic  scales.  The  same  contact 
potential, 0. 75  volts,  was  added  to  V in  computing  p.  The  cutoff  transition 
lies  roughly  above  p = 0.4.  1 A.  comes  out  of  a vertical  dive  about  at  0.04, 
and  after  cutoff  depends  not  only  on  p,  but  also  on  and  T^.  The  dependence 
on  T^  may  in  part  reflect  changes  in  gas  pressure  in  the  (sealed  off)  tube. 


! 
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Figure  4-2e  compares  the  data  in  the  four  sections  above  it.  The 
sluuB  li  ill  (I^/I j ) J ^ ill  p Is  iiearly  Iiiuepeiuleiit  of  and  decreases  as 
rises,  while  the  intercepts  oi  the  current  curves  at  p = 0.  1 increase  with 
V and  with  T . 

2k  C 

For  p below  0.  1 some  curves  of  leakage  current  (Fig.  4-3)  show 
"humps.  " Such  humps  appeared  in  the  first  published  cutoff  curves  (Fig. 

11  of  Hull's  paper  [1])  and  have  reappeared  in  subsequent  articles.  The 
hump  for  V = 1.5  in  Fig.  4-3  is  not  affected  by  T , but  the  hump  for  V = 

2k  C 2k 

6 volts  seems  to  move  toward  lower  values  of  p as  T falls.  We  have  des- 
cribed  a drag-loop  effect  of  the  probe  signal  within  one  such  hump  (Fig. 
4-5b*).  These  humps  do  not,  in  our  tube,  produce  strong  noise  or  oscil- 
lations, and  they  occur  only  well  removed  from  cutoff. 

Thus  in  summary,  if  electrons  in  orbits  of  the  emission  type  are  the 
exceptional  ones,  if  the  average  electron  suffers  many  collisions  or  a con- 
tinual exchange  of  energy  and  momentum  with  its  neighbors  during  its 
life-time  in  the  cloud,  one  would  not  expect  the  cloud  to  have  a sharp  edge 
at  8^.  Although  the  anode  current  drops  sharply  in  the  cutoff  transition, 
indeed  from  80  to  3 per  cent  of  the  Langmuir  current,  1 continues  to  flow 

2k 

for  very  low  values  of  p.  We  find  little  support  for  the  suggestion  that 
the  diede  after  cutoff  maintains  an  electron  cloud  bounded  by  a contract- 
ing radius  s^,  a radius  as  sharply  defined  as  the  cutoff  transition  is . 

d.  The  Potential  in  an  Unbounded  Cloud 

We  found  in  T.R.  185  a reasonably  sharp  and  deep  resonance  absorp- 
tion over  most  of  the  interval  0<C  p 0.  3 (Fig.  4-2b*).  We  would  have 

seen  multiple  resonances,  perhaps  smeared  into  a broad  region  of  absorp- 

2 

tion,  if  distinct  values  of  u were  characteristic  of  different  radial  intervals 

2 

'vithin  the  space-rharge  cloud.  Hence  one  value  of  u =1  - C,  for  a given 

set  of  conditions  p,  X and  T , must  characterize  at  least  roughly  the  whole 

pc  . 

interaction  space,  if  the  cloud  extends  practically  to  the  anode  after  cutoff . 

If  we  accept  a two-parameter  potential  function,  and  one  parameter  applies 
in  the  whole  radial  interval  Ks-'Cs^,  then  the  other  must  also  apply  through- 
out this  interval,  if  the  potential  ^(s)  is  continuous  in  s.  The  two  parameters 
could  be  B and  C (4.  3),  but  the  potential  may  now  be  completely  deterxnined 


r 

I 

I 


TRZOl 


-42- 


by  two  measurable  quantities,  u^  and  p/^^: 


" P <*  i\  /I  ..2vi_  j.ii  ^ 

JD  to  - i I - \ *-***•  “ ... 

p ' a ' ' 


cL 


<W8)  = - 1)  - (1  - «^)(a_  - s)3 

O S U Ab  A 


b4 

T~T 


m u r 


/A  1 ?i 


(4.13) 


K(8)  = -^[^(8,  - 1)  - (1  - U^)(8  - 6^)] 


erE 


x»  - X.  - - ,j  2 Z 

2m  u r 


(4. 14) 


aW 


1 

8s' 


Sa  - 1)  - (1 


«^)(8a  + 8^)]  = 


4m  0)  6 


(4.15) 


We  found  in  Fig.  4-2*  that  most  resonance  data  fall  near  the  line 
2 

u = P/P„ 


If  this  relation  is  added  to  (4.  12)  then  2C  + B = 0,  and  the 


next  three  quantities  are: 

*<•'  = afs- ■ KC»  = -5l5-C^-i;.  = 


At  cutoff,  p = p^,  these  formulae  arc  the  same  as  those  derived  from  Brillouin's 
model.  The  formula  for  q(s)  in  (4. 16)  was  first  used  in  (2.44*).  For  all 
values  of  p,  K(l)  = 0,  and  the  electric  field  in  (4.  16)  vanishes  at  the  cathode. 

We  regard  as  a partial  confirmation  of  our  model  of  an  unbounded  cloud  the 
fact  that  K(l)  = 0 arises  from  this  interpretation  of  resonance  data. 

In  Fig.  4-2*  points  (p,u^)  numbered  1 and  2 were  singled  out  to  represent 
the  resonance  absorption  in  our  magnetron  diode  with  a hot  cathode.  Figure 
4 '4a  shows  potential  curves  corresponding  to  these  points;  the  solid  curves 
are  based  on  (4. 13)  and  the  dashed  curves  on  Brillouin's  model.  In  each  case 
8s^  is  plotted  against  s.  Clearly  the  Brillouin  potential,  or  indeed  any 
pcteitial  based  on  a cloud  bounded  by  s^,  lies  well  above  the  corresponding 
solid  curve  for  s>  s^^. 

In  Fig.  4 -4b  the  solid  curves  are  the  standard  ones  fur  a diode  without 
magnetic  field,  the  Langmuir  and  the  logarithmic  (space-charge  free)  poten- 
tials. The  dashed  curves  reproduce  the  solid  curves  of  Fig.  4-4a,  but  the 
ordinate  has  been  renormalized  as  to  make  a comparison  of  curve 


FIG, 4-4  REPRESErrrATIVE  POTENTIAL  FUNCTIONS 
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shape  easier.  Again  the  potential  curves  based  on  (4.13)  lie  below  the  standard 
curves . 

2 

The  insert  in  Fig.  4-4b  shows  the  dot-dash  line  i - u ~ p/p  and  the 

■> 

representative  points,  1 and  2,  lying  very  close  to  this  Jine  on  the  (p,u~)-plane. 
The  dashed  line  Qq/Q  compares  the  Brillouin  model  and  the  unbounded  model 
(4.  16)  with  respect  to  total  charge  inside  the  cloud  as  a function  of  p.  Clearly 
the  cloud  extending  to  the  ancde  admits  more  charge  than  the  Brillouin  cloud 
bounded  by  Sj^: 


s 

Q(s)  = ^ r pds  = 2xh  [K(s)  - K(l)] 

4m  ut  J.  ^ 


Qb  J,  + 


(4.  17) 


Qg:  B = -2  , C = 1 , l<s<Sjj;  Q(s^):  B = -|£  , C = l<s<s^. 

*^c 

Figure  4-4c  represents  potentials  in  a diode  with  a very  large  ratio  of 

anode  to  cathode  radius,  s^.  W.  Engbert  in  1938  [13]  used  such  a aiode, 

with  filamentary  cathode  of  about  0.  2 mm  radius  and  anode  of  7 mm  radius, 

in  experimental  measurements  of  the  radial  variation  of  potential  in  a magnetron 

after  cutoff.  Kis  measurements,  the  only  ones  so  far  reported,  deserve  more 

consideration  than  they  seem  to  have  received.  The  points  © in  Fig.  4-4c 

are  repiotted  in  coordinates  (s/s  , s^/s  ) from  Engbert,  Abbild  5,  and  the 

a I / a a 

solid  curves  again  represent  the  Langmuir  and  logarithmic  potentials.  The 

dashed  curves  through  the  points  O have  much  the  same  relation  to  the  solid 

curves  as  the  dashed  curves  1,2  have  to  the  same  solid  standard  curves  in 

Fig.  4-4b.  In  particular,  the  dashed  curves  lie  well  below  the  solid  curves 

near  the  anode  and  indicate  a denser  space  charge  near  the  anode.  Engbert' s 

data  lend  no  support  to  the  models  that  confine  the  charge  cloud  inside  a 

radius  s . 

H a 

Two  figures  reproduced  from  Engbert' s paper  appear  in  Fig.  4-5. 

Abbild  2 shows  his  smooth  anode  tube  suspended  between  magnet  faces  M 
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with  cylindrical  axis  vertical.  He  could  move  a fine  tungsten  wire  radially 
through  the  interaction  space.  At  the  end  of  this  wire  he  attached  a sliver 
of  tungsten  /A).  The  potential  V^,  applied  to  the  wire  by  the  potentiometer 
circuit  shown  in  Abbild  1,  would  deflect  the  sliver  toward  anode  or  cathode. 
When  opening  and  closing  switch  T did  not  deflect  the  sliver  from  the 
vertical.  Engbert  recorded  as  the  undisturbed  potential  in  the  diode  at 
the  radial  position  r^  of  the  sliver.  Through  a microscope  trained  on  a 
slit  in  the  anode  he  observed  the  sliver  by  the  light  of  the  hot  tungsten  cathode. 
If  V was  larger  than  150  volts,  he  felt  he  could  determine  V more  ac- 
curately  than  r^.  Without  a magnetic  field  he  obtained  data  for  good  curves 
of  the  standard  Langmuir  and  logarithmic  potentials. 


After  the  magnetic  field  cut  off  the  anode  current,  he  reports,  the 

strength  of  the  field  and  the  cathode  te  mperature  T^  had,  within  widelimits, 

no  effects  On  thc  potential  distribution.  (Hi«  experimental  points  actually 

have  the  spread  indicated  in  Fig.  4-4c,  so  it  seems  likely  that  within  this 

spread  one  might  find  correlation  between  /w(r  ) and  V and  T . ) He 

P / * a c 

examined  a wnde  range  of  T^,  but  he  does  not  give  the  range  of  p he  ex- 
amined. 


Engbert' s interpretation  of  his  results  is  a little  forced.  He  approxi- 
mated his  potential  by  two  straight  lines,  meeting  at  r^.  the  radius  that  gave 
the  best  fit  to  his  data: 

V (r-r  ) 

V(r)  = 0,  r < r<r^  ; V(r)  = \ , r^<  r<r^  . (4.  18) 

a o 


In  the  first  interval  he  seems  to  believe  there  is  no  radial  motion  of  elec- 
trons; in  the  second  he  draws  electronic  orbits  like  those  proposed  by 
HiiU  [1].  Doehler  [14]  fits  Engbert's  data  with  a curve 0(rj  = V^(r/r^)*^, 
n = 0.875.  Twiss,  (5,  page  282)  referring  to  these  two  types  of  fit,  says 
"The  accuracy  of  the  measurements  cannot  be  very  high,  if  they  are  open 
to  two  quite  different  interpretations."  Perhaps  one  might  question  first 
the  accuracy  of  the  interpretations. 

For  the  present  argument  it  is  unfortunate  that  Engbert  used  a fila- 
men1:ary  cathode.  The  enlarged  corner  section  in  Fig.  4-4c  shows  how 


FIG.  4-5  ENGBERT’S  EXPERIMENT 

(HOCH  FREQUFNZ  TECHNIK  51, 
^4  -52  (1938)  ) “ 
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poorly  near  the  cathode  a quadratic  potential  function  (dot-dash  curves) 
matches  the  potential  he  measured.  Perhaps  measurements  like  his, 
but  made  on  a diode  with  a large  cathode,  would  yield  a potential  function 
to  which  one  could  fit  a quadratic  form  closely,  "fhe  reciprocalinfsrsnce, 
that  a probe  signal  resonance  in  a magnetron  diode  with  a filamentary 
cathode  woijild  be  broad,  must  be  drawn  from  Engbert's  data.  If  the  potential 
curve  has  a wide  range  of  curvature,  when  plotted  as  in  Pig.  4-4c,  then 
as  we  noted  above,  the  Individual  electronic  orbits  will  tune  to  a range  of 
frequencies.  We  have  checked  this  inference:  a pre-war  model  s^it- 
anode  magnetron  was  mounted  so  as  to  transmit  a TEM  signal  along  its 
filamentary  cathode.  In  it  the  cyclotron  resonance  absorption  was  ob- 
servable but  extremely  broad.  This  was  also  true  of  the  region  in  which 
a strong  probe  signal  tended  to  reduce  the  anode  leakage  current. 

In  summary,  then,  we  have  developed  a potential  function  for  an  un- 
bounded space-charge  cloud  in  a magnetron  diode  after  cutoff.  This  func- 
tion interprets  well  the  normal  resonance  absorption  data,  and  the  radial 
electric  field  derived  from  it  vanishes  at  the  cathode.  The  unbounded 
cloud  contains  in  fact  more  total  space  charge  than  the  saturated  clouds 

in  models  bounded  by  s„.  The  cloud  density  is  a function  of  anode  voltage 

**  1-2 
is  always  less  than  the  saturated  density,  q = -j-  (1  + » )♦  The 

potential  function  has  a family  resemblance  to  the  potential  Engbert  measured 

in  a magnetron  diode  with  much  larger  s^. 
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V 

The  Distribution  of  Electrons  in  Phase  Space 

a.  Introduction 

In  this  chapter  we  extend  somewhat  the  statistical  approach  of  Hok  [1] 
to  the  problem  of  charge  distribution  within  a magnetron  after  cutoff,  and 
consider  the  density  clist.ribution  f of  electrons  in  phase  space.  One  should, 
if  possible,  derive  i as  a solution  of  the  Boltzmann  transport  equation  (2, 
P^8®  46),  but  this  would  require  one  to  determine  the  cumulative  effects  of 
the  individual  interactions  of  electrons.  Without  attempting  that  difficult 
task,  we  present  here  in  brief  outline  some  features  of  these  interactions 
and  some  conditions  on  the  function  f . 

b.  Individual  Electronic  Interactions 

If  two  electrons  of  charge  ^in  a cloud  of  particles  approach  each  other 
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ao  closely  that  one  third  of  their  average  kinetic  energy  changes  into 
potential  energy  of  their  mutual  repulsion,  they  are  then  a distance  d 
apart: 


kT 

2 


(5.1) 


Since  the  electronic  cross  section  (or  collisions  of  this  type  is  yd / 

this  ^is  a kind  of  electron  diameter.  The  mean  free  path  h between 

2-1 

Such  collisions  is  h = (mrd  ) , if  ii  is  the  number  density  of  electrons 

in  the  cloud. 


Attempts  at  refining  this  crude  picture  of  electronic  interactions 
founder  on  the  long  range  of  the  coxilomb  force,  v/hich  tends  to  arrange 
all  the  charged  particles  in  the  cloud  in  a sort  of  lattice.  If  a section 
of  this  lattice  is  disturbed,  it  may  interact  as  a whole  with  the  electro- 
magnetic field  in  which  it  lies,  and  this  collective  behavior  may  have  a 
characteristic  angular  frequency 

Individual  particles  jostle  about  within  the  lattice  near  centers  of 
minimum  energy  in  a somewliat  random  fashion.  They  may  have  a dis- 
tribution n^  exp  |^-ed/kT]=  n^  [1  - ed/icT] , where  d(r)  depends  on  the 
distance  r_from  such  a lattice  point.  For  a similar  model  of  ions  in  a 
weakly  ionized  solution,  Debye  [3]  derived  the  equation  for  4(r): 


= 


1 d 

nr  3F 

r 


(5.2) 


The  radius  that  arises  in  (5.2)  is  called  the  Debye  screening  radius. 

It  is  the  geometric  mean  of  the  free  path  and  the  electron  radius  derived 
2 dh 

from  (5.  1):  X.^  = -j-  • The  solution  of  (5.2)  is: 


d(r)  = 


■4x6  r 

o 


(5,3) 


A weak  solution  of  ions  has  no  net  charge  and  consists  largely  of 
neutral  particles  that  act  as  a heat  bath.  Despite  stveh  differences  be- 
tween Debye's  model  and  an  electron  cloud,  BohmandPines  [4]  find  that 


TR201 


-48- 


Xjj  is  a convenient  distance  to  divide  the  interaction  range  within  the  cloud. 
Particles  which  come  closer  together  than  the  distance  they  claim, 
interact  in  a random  individual  manner,  while  particles  further  apart  than 
this  distance  recognize  only  the  collective  motions  of  the  entire  lattice  = 

The  angular  frequency  they  identify  with  the  plasma  frequency, 

2 e n 


mfe 


For  an  example,  let  us  set  e/kT  =10  volts  and  set  the  cyclotron  fre- 

10  2 -10 
quency  = 0.93  x 10  radians  per  second.  Then  we  have  X.,  = - ,65x  10 

- 3 -7 

(p  Ip)  square  meters,  h = 1.46  x 10  (p  /p)  meters  and  d = 0.89  x 10  meters. 
Moreover,  w,  — w (p/p  ),  if  from  (4.16)  we  have 

D C 


ne 

“7 — ? “ 2p 
*^c 


(1+-^) 

s 


We  have  seen  no  evidence  of  a plasma  resonance  in  our  absorption  data,  but 
one  might  expect  to  find  it  in  response  to  an  axial  electric  probe  field, 
(3.44’*').  It  would  be  quite  distinct  from  the  cyclotron  type  of  resonance, 

(1  - p/p^)  , (4.3>^). 


c.  General  Conditions  on  the  Distribution 

Individual  interactions  may  make  the  electrons  in  emission  orbits 
the  exceptional  ones  in  the  space -charge  cloud,  and  require  a statistical 
treatment.  The  statistical  model  proposed  by  Hok  [1]  is  as  extreme  in 
this  direction  as  the  single  and  double  stream,  models  arc  in  another  direc- 
tion. If  the  cathode  and  anode  and  the  ends  of  the  magnetron  diode  were 
reflecting  walls,  instead  of  source  and  sink  for  electrons.,  then  the  charge 
cloud,  he  points  out,  would  fill  the  space  between  these  walls,  and  the 
particles  in  it  would,  after  many  collisions,  approach  a thermal  eqxiili- 
brium.  The  spacial  region  of  lowest  potential  ed  (since  e,is  negative), 
wovild  be  close  to  the  anode;  and  it  would  be  the  region  of  highest  density, 
if  f were  proportional  to  the  Boltzmann  factor 

exp|-  (^^^+  ed)  (kT)  . 
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We  develop  a leas  extreme  form  £,  which  may  range  over  all  elements 
of  a ■ix-dimcnsional  phase  space  (XiytT;  r,6,z).  Here  we  use  x.^F  and  ^ - 

as  in  (2.  2).  Of  course,  ,f  is  larger  than  zero  only  within  the  cylindri- 
cal lutcraction  rcgicn  in  the  diode,  r^<;  r < r^,  and  0 <z  <h.  In  the  velocity 
subspace  (x,y,F),  some  regions  imply  orbits  which  collide  with  a physical 
boundary,  and  f should  reflect  a thin  or  one  chot  population  in  these  regions. 
Another  region  of  (x,y,z)  includes  emission  orbits  and  may  have  a denser 
population  even  though  such  orbits  carry  a small  proportion  of  the  particles 
in  the  space  charge. 

Since  tVie  leakage  current  and  the  probe  signal  data  show  for  low  a 
dependence  on  the  cathode  temperature  T^,  and  a dependence  on  that  re- 
mains after  normalizing  with  respect  to  u> , we  may  expect  f to  depend  on 
and  cn  as  well  an  on  The  first  moment  of_f,  the  charge  density  p, 

will  with  d(r)  satisfy  Poisson's  equation  (1. 16),  £ind  the  radial  velocity 
moment  of  £ will  give  the  anode  current 


(5.4) 


= 2xh  e J'JJ  xf  dx  dy  dzj^'^^ 


(5.5) 


The  domain  G of  these  integrals  is  the  domain  of  six-space  within  which ^ 
is  greater  than  zero.  W-e  have  noted  above  its  spacial  boundariee;  wc  may 
assume  now  that  f is  independent  of  the  polar  angle  0.  The  motion  of  elec- 
trons along  the  z-axis  maybe  aperiodic  and  limited  within  G only  by  col- 
lisions. We  assume,  except  close  to  the  diode  ends  in  small  regions  which 
we  ignore,  that  i is  independent  of  z.  (For  large  values  of  iT,  this  may  also 
be  untrue.) 

If  details  of  the  cumulative  interaction  processes  were  known,  one 
might  calculate  the  rate  of  transpox't  of  electrons  from  emission  orbits 
through  intermediate  orbits  to  collection  orbits.  The  over-all  rate  T/T 
is  indicated  in  Fig.  . In  general  the  electrons  are  knocked  out  of 
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emission  orbits  at  the  rate  required  to  replace  those  lost  to  the  cloud  through 
collection  orbits,  since  the  net  radial  leakage  current  is  independent  of  time 
and  of  the  radial  coordinate.  The  more  efficient  the  knocking  out  process, 
when  compared  to  the  radial  transport  mechanism  effective  throughout  the 
cloud,  the  fewer  electrons  are  required  to  complete  emissibn  orbits.  Inother 
wordb,  the  feature  cf  the  cathode  region  which  most  affects  the  interaction 
space,  namely  the  barrier  at  the  release  radius,  will  be  of  such  a height  as 
to  limit  properly  the  number  of  electrons  which  leave  the  cathode  region  and 
follow  emission  orbits  out  through  the  interaction  region. 

To  fix  the  boundaries  of  G in  the  subspace  (x,y)  we  must  know  the 
potential  function 4(s).  We  used  a qiiadratic  form  for  4>(s),  first  in  Chapter 
of  T.R.  185  to  approximate  the  potential  over  the  radial  extent  of  a single 
orbit,  and  more  recently  in  Chapter  IV  of  this  report  to  represent  the  potential 
over  the  whole  radial  interval  from  cathode  to  anode.  Clearly  this  quadratic 
form  fails  to  represent  the  cathode  region  correctly,  but  the  lack  of  multiple 
resonances  in  probe  signal  data  and  the  high  values  uf  justify  using  the 
form  in  this  chapter,  too.  These  high  values  also  indicate  tliat  electro<iic 
orbits  of  some  stability  are  typical  of  the  charge  cloud.  Hence,  although  we 
know  little  of  the  interactions  which  add  up  to  determine  f,  we  assume  they 
are  neither  so  strong  nor  so  frequent  as  to  invalidate  the  orbital  analysis  given 
in  Chapters  I*  and  II*  . 

d.  Projections  of  the  Domain  of  the  Distribution 

To  apply  these  general  conditions  to  two  specific  cases  we  choose  again 

the  typical  points  1 and  2 on  the  (p,u^) -plane.  Figure  5-1  shows  two  sketches 

of  projections  of  "safe  orbits,"  orbits  which  have  both  tuning  points 

inside  the  diode,  or  in  symbols,  *d—  velocity  coordinates 

have  been  normalized  in  terms  of  w rather  than  kT  ; 

c 


On  the  left  in  each  sketch  the  new-moon  shaped  curves  represent  extraordinary 
or  cyclotron  orbits,  and  on  the  right  appear  the  ordinary  or  magnetron  drift 
orbits.  We  met  this  distinction  in  (1.9*)  and  the  sketch  at  the  bottom  of  Fig. 


rr»r»  o A 1 Cl 
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5-1  may  make  it  clearer.  The  cyclotror.  orbits,  because  they  require 
high  kinetic  energies  (Y^'>4),  are  preauipably  rarely  occupied  although 
they  are  as  "safe"  as  the  other  type. 

Figure  5-2  shows  the  projections  of  two  sets  of  safe  ordinary  orbits 
on  the  eubspace  (?,s).  The  sets  are  based  on  potential  functions  appropriate 
to  the  same  two  points  on  the  (p,u^) -plane.  The  individual  projections  are 
plotted  from  (1.6*)  and  ( 1 . 15  * ): 


P = Vo'  Y + s ; u^  8 ^ + s (i-  - [1  - u^] ) 

g a Pj.  •* 


(5.7) 

P,  we  recall,  is  the  normalized  canonical  momentum  in  the  circumferential 
direction,  and  is  the  mean  radius  of  an  orbit  which  has  turning  points 


t®b*  ®d> 


Eliminating  P,  we  solve  (5.7)  for  Y(s),  an  expression  in  which 
the  plus  sign  represents  ordinary  orbits: 

1/2 

(5.8) 


i 

T/«>  - + 1 • ^ ^ ■ 

1 

L J 

Actually  the  dashed  lines  in  Fig.  5-2  represent  ten  ordinary  orbits,  chosen 
as  examples,  with  from  the  group  = 5,4, 3,2,  s^  =4,3, 2,1, 

The  boundaries  of  the  projection  of  G on  (T,s)  consist  of  the  vertical  lines 
s = 1,  s = s , and  the  sloping  solid  curves.  In  each  sketch,  the  upper  curve 
forms  the  envelope  of  ends  s^  of  orbits  which  graze  the  anode,  and  thelower 
curve  forms  the  envelope  of  tvrning  points  of  orbits  which  graze  the 
cathode. 

The  dot -dash  line  in  each  sketch.  Y = T(«  ),  that  cuts  each  orbit  in 

g g 

half,. is  the  locus  of  orbits  of  zero  radial  motion.  One  can  define  an  ef- 
fective potential  U by  reference  to  such  orbits,  orbits  of  minimum  energy 
among  all  orbits  having  the  same  mean  radius  s^’.  This  potential  is  plotted 
in  Fig.  5-3,  and  fOr  comparison  the  potential  X(^)  u-sed  in  Chapter  111.  The 
radial  dependence  is  given  for  the  same  typical  cases,  1 and  2.  Introducing 
the  notation  (1  - u - p/p  ) = d - 0 for  convenience,  we  have; 


U = Yg^-  8*(Sg)={^(s^-l)A4  1 + u^  (Sg-1)  - 1)  - u^8g-l)^ 


!i 


I 
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X=  - 1 - 8*(s  ) = {8 


8 


a 8-1 

■ 1)  [{-S:  . 1)  A + 1 - -2-  {-^ — 

' '■*8  P * 8 

g g 


)] 


(5.9) 


If  f were  plotted  in  a three -dimen8ional  graph  above  the  (Y,8)-plane, 

thi8  graph  probably  would  have  a ridge  or  maximum  above  the  dot-daeh 

line  Y . For  if  electrons  are  dietributed  in  8ome  statistiral  fashion  ac- 
g 

cording  to  their  radial  energies  above  the  zero  level  given  by  ll,  thie  dis- 
tribution must  tend  to  pile  up  above  Y . However,  Y is  also  the  locus 

g g 

of  the  points  s^  of  maximum  radial  motion  in  a given  orbit,  X = X . 

g g m 

Charges  in  periodic  motion  spend  more  time  near  their  turning  points 

than  near  the  orbit  center  and  this  tends  to  reduce  t?ie  density  above  Y . 

g 

Indeed,  if  /n  is  the  time-average  contribution  of  a particle  in  an  orbit 
{s^,s^)  to  the  number  density  n in  the  interval  As,  then: 


^ uu  As 


.2 , . , ^22 

s {s,Sg,jy)  = 8 X w 


(5.10) 


Formulae  (5. 8)  and  (5. 10)  suggest  that  f can  be  expressed  in  terms 
of  s and  j , instead  of  X and  Y.  This  cannot  be  carried  out  entirely  while 
the  form  for  the  unbalanced  radial  transport  current  is  unknown,  but  for 
the  residual  part  of  the  distribution  ure  suggest  the  following  semi- 
Boltzmann  form: 


g (b  ) 
»o'  g' 


exp{-b(s^)(J^+j^)} 




- (s* 


dj  ds  dj  ds* 
•^r  g •'z 


(5.11) 


Here  g^t^g)  is  partly  a normalization  factor.  In  the  Boltzmann  exponential, 

b(s^)  represents  a reciprocal  temperature  or  characteristic  energy  and 

depends  no  doubt  on  T and  V . In  the  same  exponential  term,  j = 

-1  c a r 


ttX 


m 


A 


Z -W  CA  * • 2 

is  the  radial  action  and  j is  a renormalization  of  (z)  . 

±±i^  * rr 

The  square  root  in  the  denominator  represents  J/s  . The  variable  s*  in  it 

is  to  be  integrated  over  the  orbit  interval  between  the  limits: 

. IVi  + ^+ 

g 


“b'^d  = “ 


ne  uruit  inw£V«u.  u:si.we 


(5.12) 


Equations  (4.  16)  and  (5.11)  provide  a hypothesis  which  can  be  tested 
further  in  several  ways.  For  instance,  one  might  repeat  the  temperature 


- \i  V- 


C' 


\\ 

\\\ 


PROJCCTION  OP  SAPE  ORBiTS  ON  TMi:  (»,Y)  PLANE 

2r^  9 

Y«- — M'/fcT 
'c- 

«EN\CLOPP  OP  NNER  TURNING  f’OiNTS  OP  ORBiTS 
^ WHICH  REACH  $ • $_ 


line  op  circular 

ORlPT  ORBITS  AND 


excludes  \ ^ 

v*0  BECAUSE  ^ ^ 

potential  min. 

LIES  IN  FRONT  CP 

cathode 

ENVaOPE  OF  OUTER  . 
TURNING  POINTS  OP 
ORBITS  WHICH  START 
ST  %•  I 


P *OJ  . u‘.Q7 


^ENVELOPE  OF  IF 


ENVELOPE  Of  »„  IF  «j.l-- 


Z>  P»0lJ,u*.a3 


FIG,5-?F'ROJtCTION5  Of  OAFt  CnSITS  ON  THE  ts,’r)  PLANE 


FiG  5-3  EFFECTIVE  POTENTIALS 


TR201 


-5  3" 


meaBixrements  (Linder  and  Wigdortschik)  and  electrostatic  probe  experi- 
ments (Engbert)  mentioned  in  Chapter  IV,  and  even  Wasserman's  [5]  stady 
of  the  magnetic  field  of  the  charge  cloud,  with  refinements  suggested  by 
this  analysis.  The  interference  of  the  TEM  probe  signal  with  the  trans- 
port process  and  the  observed  reduction  in  anode  current  at  absorption 
resonance  reqviires  further  study.  Finally  recent  work  in  France  on  the 
space -charge  amplification  of  noise  in  a linear  magnetron  ("dlocotron") 
bears  directly  on  the  distribution  of  space  charge  after  cutoff^  6,7] . 
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